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Introduction 



Recent theoretical and experimental efforts in quantum optics and quantum information have been focused 
on the engineering of highly nonclassical, non-Gaussian states of the radiation field IT], in order to achieve 
either enhanced properties of entanglement or other desirable nonclassical features lU |2l [3] IH |5| . It has 
been shown that at fixed covariance matrix, some of these properties, including entanglement and distillable 
secret key rate, are minimized by Gaussian states JS). In the last two decades increasingly sophisticated 
schemes for the generation of non-Gaussian states have been proposed, based on delocalized photon addition 
or subtraction iTl [8l l9l [TOl fTTIl : or on strong cross-Kerr interactions lfT2l . Some of the photon addition and 
subtraction schemes have been experimentally implemented to engineer non-Gaussian photon-added and 
photon-subtracted states starting from Gaussian coherent or squeezed inputs ifTSl [T4l [TSl . 

Through the photon subtraction of a single delocaUzed photon from a two-mode entangled initially Gaus- 
sian state, it has been possible to realize an state of enhanced entanglement with negative two-mode Wigner 
function lfT6l . Remarkably, the photon-addition/subtraction operations, performed on thermal light fields have 
led to the demonstration of the commutation relation rules for quadrature operators; one of the constitutive 
relations of quantum mechanics ifTTl . Moreover, very recently, a protocol has been experimentally realized 
that allows the generation of arbitrarily large squeezed Schrodingercat states, using homodyne detection and 
photon number states as resources ifTsl . This class of optical cat states is of particular importance because it 
is strongly resilient against decoherence lfT9l . 

Progresses in the theoretical characterization and the experimental production of non-Gaussian states 
are being paralleled by the increasing attention on the role and uses of non-Gaussian entangled resources in 
quantum information and quantum computation with continuous-variable systems ll20l . Concerning quantum 
teleportation with continuous variables (CV), the success probabihty of teleportation can be greatly increased 
by using entangled non-Gaussian resources ll3l 1211 l22l |23] l24l . In refs. 13] |2T1 |22l |23l conditional measure- 
ments, inducing "degaussification" through photon-subtraction, are exploited to improve the efficiency of CV 
teleportation protocols. Most of these investigations have used the transfer operator ll25l l26l l27l formalism 
and Fock basis representations of the degaussified resources rather than phase-space representations as used in 
the original CV protocol ll20l . Phase-space and conjugate phase-space representations of operators constitute 
an unifying language for the description of quantum optical states and processes; including the composition 
of said processes into protocols such as CV quantum teleportation. 

Moreover, non-Gaussian cloning of coherent states has been shown to be optimal with respect to the 
single-clone fideUty 1281 . Determining the performance of non-Gaussian entangled resources in CV quan- 
tum communication protocols can prove to be useful in a number of concrete applications ranging from 
hybrid quantum computation ll29l to cat-state logic l30l and in all quantum computation schemes based on 
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communication that integrate together qubit degrees of freedom for computation with quantum continuous 
variables for communication and interaction 1311 . 

In the same fashion, but following a more general approach, in ref. 1241 the implementation of simulta- 
neous phase matched multiphoton processes and conditional measurements are used to introduce a general 
class of two-mode non-Gaussian entangled states, in the form of squeezed Bell-like states endowed with a 
free parameter The optimization on this free parameter allows a remarkable increase of the teleportation 
fideUty for various classes of input states. 

In the present work, we propose a general formalism for the study of Quantum Teleportation in a CV 
setting, based on the Wigner's characteristic function description as a conjugate phase-space representation 
and on the Weyl's correspondence ll32l of unitary evolutions and measurements to coordinate transformations 
and partial traces over the characteristic function. We show that the entire quantum teleportation protocol as 
formulated in 1331 for any combination of entangled resource state and of input state can be represented by 
adequate transformations and integrations of the characteristic function of the joint input-resource physical 
system on which the protocol is performed. We formulate a number of complications to the original protocol, 
intended to simulate imperfections in the homodyne detection apparatus and the presence of environmental 
"noise" in the preparation setup for the resource state l34l[35l . 

With the theoretical tools above described, we investigate systematically the performance of different 
classes of entangled two-mode non-Gaussian states used as resources for continuous-variable quantum tele- 
portation. In our approach, the entangled resources are taken to be non-Gaussian ab initio, and their properties 
are characterized by the interplay between CV squeezing and discrete, single-photon pumping. Our first 
aim is to determine the actual properties of non-Gaussian resources that are needed to assure improved 
performance compared to the Gaussian case. At the same time, we carry out a comparative analysis between 
the different non-Gaussian cases in order to single out those properties that are most relevant to successful 
teleportation. Finally, we wish to understand the role of adjustable free parameters other than squeezing, in 
order to sculpture resources to achieve optimized performances within the set of non-Gaussian resources. 
With this objective in mind, we propose the squeezed Bell-like states l24l with a single, free superposition 
parameter which determines non-Gaussianity. The squeezed Bell-like states are formulated to include as 
special cases all the other Gaussian and non-Gaussian resources evaluated in this work. Then, we optimize 
the teleportation fidelity using squeezed Bell-Uke resources with respect to the superposition parameter We 
show that maximal non-Gaussian improvement of teleportation success depends on the nontrivial relations 
between enhanced entanglement, suitably measured level of non-Gaussianity, and the presence of a proper 
Gaussian squeezed vacuum contribution in the non-Gaussian resources for large values of the squeezing 
parameter (squeezed vacuum affinity). 

The squeezed Bell-like state can be parameterized as a first-order (or two-photon) truncation of the 
squeezed vacuum states. A generalization of the squeezed Bell-like state can be easily constructed by 
allowing four-photon terms in the (prior to squeezing) superposition making up the Bell state, thus con- 
structing a more general superposition of Fock states l37ll . Further optimization of teleportation fidelity is 
made possible by the addition of a new dimension to the Hilbert space to be explored and the addition of a 
new free parameter to the set of parameters over which optimization is carried out. However, optimization 
reduces these squeezed superpositions of Fock states to second-order truncations on squeezed vacuum states. 
An avenue of exploration of non-Gaussian resources other than superpositions of few-photon Fock states 
lies in the formulation of resources combining two-mode squeezing and the entangled superpositions present 
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in two-mode Schrodinger Cat states. Such two-mode squeezed cat-like states can likewise be used in CV 
teleportation and optimized with respect to the free parameter given by the phase-space distance between the 
terms of the superposition. 

Furthermore, we investigate the effects of the presence of thermal noise on the performance of two- 
mode non-Gaussian states used as resources for continuous- variable quantum teleportation ||36l [37l . We will 
consider the non-Gaussian resources obtained by superimposing the classes of squeezed Bell-like states and 
squeezed cat-like states over two-mode thermal states 1381 . Due to the thermal contribution, the state so 
obtained is mixed and its correlation properties are modified and deteriorated by the presence of thermal 
photons. We limit the discussion to the situation of ideal teleportation protocol, i.e. ideal Bell measurements 
and decoherence-free propagation through space of radiation states. Detailed analysis in the instance of the 
most general realistic situation, including various sources of noise, will be discussed elsewhere. 

The thesis is organized as follows. In chapter[T|we discuss the basics of CV systems (see section fTTTT i. i.e. 
quantum optics including, most importantly the groups of operations such as displacement and squeezing and 
non-unitary "operations" such as homodyne measurement (see section [L2] |. We introduce correspondence 
principles based on the aforementioned groups of operators and phase-space representations of states of the 
radiation field, as well as the characteristic functions of some Gaussian states (see section [T3T l. Finally, we 
discuss quantum teleportation as an universal procedure and teleportation fidelity as a measure of teleportation 
success with a view to the formulation of quantum teleportation in a characteristic functions' language (see 
section [L4l i. 

In chapter|2]we derive, using the Wigner's function (see section [OT i and Wigner's characteristic function 
language (see section [ZZl i, the expression for the output state of CV teleportation for the ideal setup and for 
several interesting modifications of the teleportation protocol. Such as the performance of the "homodyne" 
projective measurement over a mixed state (see section I2.3l l: a resource state prepared and superimposed 
over thermal vacuum states (see section |24] i and "realistic" homodyne detection (see section [23] ). whereas 
fictitious beam-splitters mix the modes to be measured with external fields and cause a loss of intensity). We 
analyze and compare the effect of the complications and modifications introduced in each case. Finally, we 
analyze the teleportation fidelity expression we have chosen |f39l |40| for the teleportation outputs we have 
derived in the characteristic function formalism(see section l2!6l l. 

In chapter [3] we study the use of some non-Gaussian states as resources for an ideal CV teleportation 
protocol. We introduce and describe relevant instances of two-mode entangled non-Gaussian resources, 
including squeezed number states and typical degaussified states currently considered in the literature, such 
as photon-added squeezed and photon-subtracted squeezed states (see section [TTI ). We compare the relative 
performances of non-Gaussian and Gaussian resources in the CV teleportation protocol for different (single- 
mode) input states, Gaussian and non-Gaussian, including coherent and squeezed states, number states, 
photon-added coherent states, and squeezed number states (see section 13.21 ). We introduce the squeezed 
Bell-like states as a generalization including all of the former non-Gaussian resources, as well as the Gaus- 
sian two-mode vacuum and squeezed vacuum (twin-beam) as special cases (section 13.31 ), and consider the 
optimization of non-Gaussian performance in CV teleportation with respect to the extra angular parameter 
of squeezed Bell-like states, and show that maximal teleportation fidelity is achieved in every case using 
a form of squeezed Bell-like resource tailored to the input that differs both from squeezed number and 
photon-added/subtracted squeezed states. We identify some properties that determine the maximization of 
the teleportation fidelity (see section[331i using non-Gaussian resources; finding that optimized non-Gaussian 
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resources are those that come nearest to the simultaneous maximization of three distinct properties: the 
content of entanglement, the amount of (properly quantified) non-Gaussianity, and the degree of "vacuum 
affinity", i.e. the maximum, over all values of the squeezing parameter, of the overlap between a non- 
Gaussian resource and the Gaussian twin-beam. Schemes for the experimental production of optimized 
squeezed Bell-like resources are proposed and illustrated (see section [375] l. 

In chapter |4] we introduce a higher-order generalization of the squeezed Bell-like states: squeezed su- 
perpositions of Fock (SSSF) states and a new class of resources, the squeezed cat-like states. We define, 
study and optimize the new class of SSSF states (in section 14.11 ). We show that all the squeezed Bell- 
like states and the optimal SSSF states can be regarded as "truncations" on Gaussian states. Higher order 
"truncations", by bestowing an extra dimension to the Hilbert space on which optimization is performed, 
further improve fidelity, over the already optimized fidelity of squeezed Bell-like resources. We introduce the 
cat-like resources, two-mode squeezed superpositions of coherent states (see section |4!2] |. We optimize the 
cat-like states for fidelity of teleportation: We find them to be non-Gaussian resources with a teleportation 
performance inferior to that of the optimized squeezed Bell-like state; but nevertheless superior to that of 
the Gaussian states for equal squeezing. In the two sections of this chapter, we perform an analysis of the 
entanglement, non-Gaussianity and Gaussian affinity for all the resources, similar to the analysis of the same 
properties performed in chapter[3] 

chapter|5]refers to the teleportation protocol using the general class of squeezed Bell-like states (and thus 
all the non-Gaussian resources introduced in chapter [3]), together with the cat-like resources (introduced in 
section Wl% in the presence of noise; for the teleportation of coherent state inputs. The resource has been 
prepared or propagated, i.e. superimposed in a noisy environment made up of thermal states, resulting in 
a mixed-state resource ll38l . First, we compare the performance of the mixed squeezed Bell-like states and 
mixed squeezed cat-like states when optimized for maximum fidelity and a similarly mixed Gaussian resource 
(see section lSTT l. Thus we study the simplest instances of teleportation using mixed non-Gaussian resources. 
We compare the robustness of the entanglement of squeezed Bell-like states, squeezed cat-like states and two- 
mode Gaussian states under noisy conditions (see section fOl ). Firstly, considering the violation of a sufficient 
inseparability criterion 11411 l42ll for mixed Bell-like and mixed Gaussian states at given levels of noise. Lastly, 
by considering the noise-induced arrival of the teleportation fidelity at the classical teleportation threshold 
for coherent state inputs as a practical criterion for the disappearance of the entanglement of teleportation 
when the resource is noisy: for squeezed Bell-like, squeezed cat-like and Gaussian resource states alike. 

In chapter |6l we present the conclusions of our work and discuss the possibility of extending the char- 
acteristic functions' formalism to more general teleportation setups; of optimizing the CV protocol itself for 
non-Gaussian resources and inputs; and of considering other non-Gaussian resources for the analysis and 
optimization of teleportation performance. 
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Initiation 



In this chapter we will go through some basic and previous concepts that are necessary to the understanding 
of the main body of this work. We will also establish some conventions and definitions that will hold for the 
next chapters. 

In section fTTI we review the basic concepts of the Continuous Variables (CV) representation of quantum 
states of the radiation field, which are the subject matter of Quantum Optics. In section [TT2I we recall linear 
transformations on Continuous Variables and the procedure of homodyne detection together with some quan- 
tum states associated with these operations. With the purpose of making clear the correspondence principle 
between density operators of quantum states and phase-space (displacement operator) representations of such; 
and of establishing a correspondence between the transformations ( ideally associated with experimental 
procedures) on density operators and coordinate transformations on phase-space representations. 

In section [T31 we review the Wigner function and the Wigner characteristic function; phase-space (and 
conjugate phase-space) representations that associate an square-integrable operator (particularly the density 
operators) with an analytic function of a complex variable that functions as a pseudo-phase-space coordi- 
nate. An special emphasis will be made on the conjugate phase-space representation given by the Fourier 
Transforms of phase-space functions; the characteristic functions corresponding to observables and quantum 
states. 

We explain in section [T~4l the basic concepts of entanglement, maximally entangled state, and universal 
quantum teleportation; for physical systems with state vectors belonging to an arbitrary Hilbert space. Lastly 
we analyze, briefly, the definition of teleportation fidelity we have chosen for this work. 



1.1 Quantum optics and continuous variables 

The quantized electromagnetic field has a field operator for the photon particle; which is the electric field (or 
the magnetic field, depending on the choice of phase). For a single frequency ujk and a single polarization 
component the electric field reads 
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= E+{f,t) + E^{r,t) (1.1) 



The correlation functions of the field are the mean values of normally ordered products of E'^ and 
E^ for the appropriate times ti,t2, ■ ■ ■ and positions ri, r2, • ■ •■ For example, the second-order correlation 
function for the field (frequency oj^.) is given by 



{E-{r^,t')E+{r,t)) (1.2) 

where, for = r and t = t',we have the field intensity, or average number of photons with energy hujk for 
position r and time t. 

The Hamiltonian of the radiation field, which for a classical field is given by 



H = 2-' j dr [\E\' + \B\'j (1.3) 
becomes, for the quantized field and in terms of the annihilation and creation operators; 

Hk=Yl hojkialdk + (1.4) 

k 

where the sum is over all the light frequencies allowed by the boundary conditions established beforehand 
for the quantization of the field. 

The creation and annihilation operators for the field are those of a harmonic oscillator with bosonic 
excitations. These have constitutive relations given by their commutators; 

[dk, cil,] = 4,fc' 

\dk , dk'] = [dl , dl,] =0 (1.5) 

We will limit our review in this work to one frequency only; barring the rigorous study of non-linear 
quantum optical phenomena such as down-conversion, the generalization to multiple frequencies is straight- 
forward. The Hamiltonian of eq. ilAi . limited to one frequency uj, is linear in the number operator of the 
harmonic oscillator fi = d'^d. 

The eigenstates of the number operator have a definite number of photons ll44ll . They are called the 
number, or Fock states In ); 



a\m) ^ \Jm\m ) 

|m ) = \/m + 1 1 771 ) 
n\m) = m I m ) 

(n I m) =Sn,m (1-6) 



' In the sense that aj. is always to the left of aj; 
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for m = 0, 1, 2, . . . 

We can define the position and momentum operators of the harmonic quantum oscillator; treating it like 
a particle in a quadratic potential, with unit mass. Define the annihilation and creation operators as 



1 , 

(uxk + ip) (1.7) 



a) = —j= {ujx — ip) (1-8) 
H = ^{p^ +Lo^x^) (1.9) 

The position and momentum operators of the harmonic oscillator of unit mass will correspond to the 
quadrature operators of the radiation field. The eigenvalues associated with these Hermitian observables are 
real numbers; the continuous values of position and momentum, and the reason for the Continuous Variables 
name given to systems thus describable. 

With the commutation relations of eq. ( 11.5b we can calculate the commutation relations for the quadrature 
operators, 

[x,p]^ih6 (1.10) 

that define, in turn, the Heisenberg uncertainty relation obeyed by these operators. Given Ax = x — (x) and 

Apk = p — {p) we have; 

{Ax') {Ap') > \\{[x,p])\' = ^ (1.11) 

For simplicity, we choose the system of units whereby h = 1/2 and uj = 1. In this way the quadratures 
become dimensionless; 



X ^ - {a + a)) = Rc[a] 

p ^^.{a-a)) = Im[a] (1.12) 

In this manner; a = x + ip and the physical quantities associated with observables x and p can be 
made to correspond with a complex phase-space "coordinate" a, for the representation of CV states and 
operators 1321 . The Hermitian, dimensionless position and momentum quadratures would correspond, re- 
spectively, to the real and imaginary parts of the coordinate a = x + ip. For example, the average values of 
(x) and {p) correspond to the average phase-space "coordinate" of the quantum state. 

However, the Heisenberg uncertainty relation (eq. (11. lit ) precludes the joint knowledge of the values of 
momentum and position; x and p with arbitrary precision for a given quantum state. This makes it impossible 
to define a as a genuine phase-space coordinate with definite values; this is just not allowed by quantum 
mechanics. In an analogous manner, having a associated as a physical quantity to an observable operator 
having an orthonormal basis of eigenstates is impossible; the corresponding operator a is not Hermitian. The 
classical radiation field is not subject to such a fundamental constraint on the precision of the joint knowledge 
of its phase and amplitude; its phase-space localization. 

We can however, define the eigenstates of the observable operators x and p. These are the position and 
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momentum eigenstates \x ) and \p). The position eigenstate \x' ), for instance, will have a position x' which 
is unambiguously defined; {x'\ Ax^ \x') = 0. For this state, the Heisenberg uncertainty relation of eq. ( 11.1 11 1 
forbids any precision in the knowledge of momentum for such a state; for {x' \ Ap^ \x') = oo. In an analogous 
manner, the momentum eigenstates are of known momentum and undefined position. 

The position and momentum eigenstates form an orthogonal basis for the representation of quantum 
states fi4l . The wave functions of a quantum state of vector can thus be defined as the projections 
ipix) = {x\'>p) and '>p{p) = {p\ip)- Given that {x\p) ~ -k^^I"^ e?^^^ , the wave function in the momentum 
representation is the Fourier Transform of the wave function in the position representation; 



The relation between wave functions in position and momentum means that a wave function narrow 
in the position representation (small ( Ax^ )) will be wide in the momentum representation; following the 
Heisenberg uncertainty relation of eq. dl.l lb . For the position eigenstate itself, the position representation 
wavefunction will be of the form {x\x') = 6{x — x'). 

The position and momentum eigenstates, useful for representing quantum states are, however, physically 
unfeasible. The average number of photons of these states is infinite; an infinite amount of energy would be 
needed for their preparation. This can be easily seen, as ( Afi^) = oo where An = n — (n) for any position 
or momentum eigenstate. The wave functions of the position and momentum eigenstates are therefore not 
square integrable; as the energy of a wave is equal to the integral of the square of its modulus. 

However, physically feasible approximates of the position and momentum eigenstates exist, they are the 
squeezed states referred to in subsection |1.2.4| 

Define nonlocal quadratures for a two-mode state, that are linear combinations of the quadratures of 
two (one-mode) states; we have as the eigenstates of such nonlocal quadratures the well-known Einstein- 
Podolsky-Rossen (EPR) states |45 |. 

In section n~.2.2| we will see how the EPR state comes about from the mixing of a position and momentum 
eigenstate by means of a beam-splitter transformation. For the same reasons given for the position and 
momentum eigenstates, the (EPR) states have infinite energy and are physically unfeasible. 

1.2 The toy box: beam-splitting, squeezing, displacement and homo- 
dyne detection 

We will review in this section the unitary transformations, acting on one and two modes of the radiation field, 
that constitute a basic toolbox of CV transformations and a basis for the representation of density matrices|f|of 
CV states; together with the bases of pure quantum states directly associated with these transformations. We 
will also describe the projective measurement of a quadrature of the radiation field by homodyne detection. 

1.2.1 The displacement operator and the coherent states 

The displacement operators 13811461 and the coherent states 1381 associated with them constitute the basis of 
the conjugate phase-space representations we will use in this work. 

-and any bounded operator F having a finite Hilbert-Schmidt norm Tr (hat F'^ F) 




(1.13) 
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Define the displacement operator; 

i5(a) = e(""^-"*^) (1.14) 

where a is a complex number 

The displacement operators form an unitary, orthogonal group of transformations under the operator 
multiplication map and the form trace of a product of operators. Given operators A and B with general 
ordering identities ||47ll 

e^Be-^^B + [A,B] + ^[A, [A,B]] + ^[A, [A, [A, §]]] + ■■■ (1.15) 
e^e^ ^e^+^+i/^I-^.s] [[A,B],A]^[[A,B],B]^0 (1.16) 

and the commutation rules for annihilation and creation operators of eq. ( 11.5b . it can be shown that | 



D^a) ^ D-\a) = D{-a) (1.17) 
D{a)D{p) = e ^ ("'''-"''''^(a + (3) (1.18) 
TY{D{a)D-^ {p)) ^TT5'^^\a~p) (1.19) 



The states most obviously associated with the displacement operator are the coherent states P8| produced 
by the displacement transformation of (multiplication by the displacement operator) the vacuum state of the 
quantum harmonic oscillator |0). The coherent state \a) is an eigenstate of the annihilation operator a with 
complex eigenvalue a. The coherent states form a non-orthogonal, over-complete basis of representation for 
one-mode states of radiation; 



l«>= D{a) |0) 
a\a) = a \a) 

oc 

\a) = e^'""' ^ a'^n!^"-^/^^ \n) where { |n )} is the Fock states basis. 

n 

cfa \a) (a\ = i 



{P\a) 



|a|V2-|/3|V2+/3*c 



(1.20) 



The averages of the dimensionless position and momentum observables are (x) = Rc[a] and {p) = 
lm[a] for coherent state | a) , with the minimum uncertainty allowed by the Heisenberg uncertainty relations 
(eq. dl.l lb ): ( Ax^ ) = ( Ap^ ) = 1/4. Furthermore, the wave functions of the coherent state on the 
momentum and position representations are Gaussian; therefore completely defined by the average values 
above. 

Coherent states are therefore the closest approximation allowed by quantum mechanics to definite local- 
ization in phase-space; they are defined by their average "phase-space coordinate" a. In the course of this 
work this "coordinate" will be given by a = x + ip; bearing in mind that a genuine phase-space coordinate 
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with definite values has no reaUzation in quantum mechanics. 

The coherent states are the closest approximation to a classical radiation field of known phase and 
amplitude among the quantum states of radiation: Because of the near-definite localization in phase-space, 
which uniquely identifies each coherent state as it would a coherent, classical field of one mode or a point 
particle; and because of their optical coherence properties 1381 . 

The transformation effected by the displacement operator on the one-mode annihilation operator is straight- 
forward to derive, bearing in mind the operator ordering identities of eqs. ( 11.151 ). ( |1.16t and the commutation 
rules of eq. (11.51 1 

D{a)aD''{a) ^ a + a (1.21) 
Lastly, using the ordering identities of eqs. (11.15b . ( 11.16b . together with the fundamental identity 

5(2)(a) =^-2yd2^e"«*-"*« (1.22) 

and the completeness properties of the coherent state basis (see eq. ( 11. 2. lb ) it can be shown that the displace- 
ment operators [46| are a complete, orthogonal basis for the representation of arbitrary bounded operators. 
Let A be bounded; such that it's Hilbert-Schmidt norm || A|| — Tr{A''A) is finite. There exists an one-to-one 
correspondence between the bounded operator A and the square-integrable form Tr{AD{^)) such that 

A= TT-^ J d^^TiiAD{0)D-\0 (1.23) 

The set of Displacement Operators {D{a), VaeC} form an orthogonal group under multiplication and 
under the trace of the product of two operators, and constitute a complete basis of representation for opera- 
tors acting on quantum states, while coherent states constitute an over-complete basis of representation for 
quantum states. These properties of both Displacement Operators and coherent states are the mathematical 
foundation for the correspondence ll32l of density operators of quantum states (and thus quantum states) onto 
functions of phase-space "coordinates" such as the Wigner Function B6l l48l and the Wigner characteristic 
function. 



1.2.2 The beam-splitter transformation and nonlocal states 

We describe below the transformation effected by an idealized linear optical device; a lossless, phase-free 
beam splitter on the two modes entering its ports. The beam-splitter transformation on the modes' annihilation 
operators is unitary; it preserves the commutation relations (see eq. ( 11.5b and overall photon number between 
incoming and outcoming modes. 

Given two incoming modes represented by their annihilation operators in the Heisenberg interaction pic- 
ture; Ga and Gb, and the two outcoming modes' operators illustrated in fig. (11.2.2b ; the transformation 
effected by a lossless, phase-free beam-splitter is given by ||2l l49l 

a„ ^ Bab{&) ab Bl^iQ) 
Bab{Q) = e^i^'>'^-'-^^-^'') (1.24) 
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a. 








a. 



Figure 1.1: Geometric configuration of an ideal beam-splitter: The incoming modes (dashed line) and 
(dotted line) are "mixed" by means of a partially reflecting surface of transmission coefficient cos(8). The 
outcoming modes are au and Oi,; each consists in a linear combination of the reflected fraction of one mode 
and the transmitted fraction of the other mode. 



where the transmittance and reflectance coefficients of the beam-splitter apparatus are, respectively, cos^(8) 
and sin^(0). Given that I3ab is an unitary transformation the commutation relations of a„ and with their 
Hermitian conjugates and between their associated quadratures will be those of Oa and ai, and associated 
quadratures. The overall number of photons is conserved, as hu + fiu ^ fia + fib- 
Let us recall eq. (11.15b : it is straightforward to re-state eq. ( 11.24b in the matrix form 



cos(0) — sin(6)\ /da 
sin(e) cos(e) / Idf, 



(1.25) 



Given that ciu.i 



ipu.v, the outcoming modes' quadratures will be linear combinations of the 



incoming modes' quadratures with an analogous relationship to that of eq. ( |1.25b . 

An operator that is an analytic function of the modes' operators, such as F{aa,a,i,; a\,al) will be 
transformed by eq. (11. 25b onto 

F'(d„, d„; al^al) Bab{Q) F b\^{Q) = F{aa{a.u, Uy), db(a„, di,); d^ (d^ , a]^), d|^(dt , d^ )) (1.26) 
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where aa^biau, ^v) denotes the inverse transform of eq. (11.25b . 

Most importantly, displacement operators for two separate modes Da{cta) Di,{ab) will be transformed to 

Dai^aiau, ay)) Db{ab{au, ay))) = Dy{ay). Where 

au\ ^ ( cos(e) - sin(e)\ / ac\ ^7) 
ayj I sin(9) cos(0) / I at / 

A separable two-mode state of radiation with a density matrix pa ® Pb entering the beam-splitter will be 
have, after the beam-splitter, a density matrix depending on modes' operators a„ , ay and their Hermitian 
conjugates. The resulting state will be usually JJ] entangled, as the density matrix will not be factorized into 
two separate density matrices for the outcoming modes u and v. 

Consider for simplicity's sake the separate wave function for a pure state il'a{xa)''Pb{xb)', after the beam- 
splitter transformation it becomes (according to eq. ( |1.25l l) equal to ipa{xa{xu, Xy), Xb{xu, Xy)). The limit 
case for the entangled states achievable by beam-splitter interaction in quantum optics illustrates the point 
nicely; assume ipa{xa) = e^'^"^ , a position eigenstate, and ipbixb) = 5{xb — x') a momentum eigenstate. 
After the beam-splitter transformation the joint wave function is 

(5((cos(e)x„ - sin(e)x-„) - x') e2''(™^(e) + ^")p' (1.28) 

which, for 8 = 7r/4 and for the modes u and v is the wave function of the maximally entangled state in a CV 
setting, the EPR state 1451 . That this state is maximally entangled can be seen easily; it is a joint eigenstate 
Hof the continuous, nonlocal quadratures 2^^/^(x„ — Xy) and 2~^/^(p„ + Py), with eigenvalues x' andp', 
respectively. While the measurements effected on one of the local quadratures, say Xu, will yield a random 

result x^u^\ with a constant probability for all the values of Xu', this same measurement will fix the value of 

^ „ 91/2^,' _ ^(™) 

It has been shown that the EPR state is physically unfeasible, unless for an infinitesimal normaUzation 
constant, the wavefunction of eq. (11.281 ) is also not square-integrable. 



1.2.3 Homodyne detection 

The procedure for homodyne detection of quadratures of the radiation field is the basic detection scheme 
of CV quantum information protocols; such as quantum teleportation ll33l l50l ISTl l52l |53]| and quantum 
tomography 1541 |55l . We will describe the experimental procedure for homodyne measurement and the 
simple projective measurement that (ideally) projects the mode thus detected onto an eigenstate of the 
measured quadrature. 

The experimental scheme for balanced homodyne detection l56l is illustrated in fig. ( I1.2.3I I. 

Mode ai is mixed with a reference mode 0.2 by means of a symmetric beam-splitter Mode 0.2 is in a 
coherent state of a very high average photon number (light intensity), approximating a classical coherent 
source of light; it's behavior can therefore be described approximately by it's complex amplitude, thus 02 ~ 



^the associated observable operators commute: 2 — Xu,Pu + Pv] = 
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Figure 1.2: Idealized Homodyne Detection Scheme: Mode ai, to be measured, is mixed by means of a 
symmetric beam-splitter with mode 0,2, prepared in a coherent state \a) of high intensity jap. Photo- 
detection is performed on modes a[ and Oj. The photocurrents ii and Z2, proportional to the measured 
intensities, are subtracted to obtain a difference photocurrent that is proportional to the quadrature's value. 

a = jaje**. The outcoming modes of the beam-splitter read 

a'2=2-^/^{\a\e"^ - ai) (1.29) 

The intensity of the two outcoming modes is measured by photo-detectors (see fig. (11.2.3b ). resulting in 
the photo-currents ii = g{a-^a[) and i2 = 5(02^02)' which are proportional (with a known gain factor g) to 
the number of photons detected. Subtracting the photo-currents, we obtain 

Si= ii-i2-- g|a|(e-'*ai +e**a];) (1.30) 

The observable measured in eq. ( II. 301 is equal to xi for $ = and is equal to pi for $ — it/2. We can 
measure, controlling the phase of the reference coherent state a generalized quadrature of mode 1; 

xl*)= 2-i(e-**ai +e'*al) (1.31) 

which would be identical to ii, for a phase shifted mode 1 where ai aie~**. The conjugate quadrature 
to xi, satisfying the commutation relations in eq. (I1.5l is simply 2;*^'^^^. 

Ideally, the measurement of the quadrature ij*'' on an arbitrary state pi of mode 1 will obtain a random 
measurement result a;'*-' and "collapse" the state into the pure eigenstate |a;'*^)i ofx^*"*. With a probability 
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equal to 



Tr(pi|x(*))i(zW|i) 



(1.32) 



thus, the system after the measurement can be said to be in a mixed state Q The pure components of the 
mixture being the Ix'*-*)! eigenstates, with probability given by eq. ( 11.32b . 

To illustrate the concept of a projective measurement further, and particularly for two-mode states, let us 
consider a state with a density matrix pi^a. The state of mode 3 after a projective measurement on mode 1 
obtaining a result a;'-*-' can be represented by a partial trace of the form 



where P(a;(*^) is the probability density for the measurement result x^*' when the quadrature x^" is mea- 
sured and y53(a;(*' ) is the state of mode 3 after measurement. The operator of eq. ( |1.33| is not a normalized, 
proper density operator, because the "projection" is not an unitary operation. It is easy to see also that the 
state of mode 3 will not be affected by the measurement process on mode 1 if the pi.s state is separable into 
two density matrices; pi, 3 = pi (g) p3. If pi, 3 is an entangled, inseparable state; the state of mode 3 after 
the measurement depends on the random outcome x'*' of the measurement performed on mode 1, with a 
probabihty P(a::(*') for the state p3(x''*-'). Such a state of "classical ignorance" produced by measurement 
is a mixture of pure states. Therefore, we choose a properly normalized state for mode 3, accounting for the 
random outcome of a projective measurement; 



1.2.4 The one and two-mode squeezing operator and the squeezed states 

Producing physically feasible states approximating as an asymptotic limit the position (or momentum) eigen- 
states, that can then be entangled by a beam-splitter requires interactions of a finite energy that cause these 
states to have a narrower position (or momentum) wave function and a wider momentum (or position) wave 
function. Though any state that has a different variance for position and momentum can be thought of as 
squeezed, we will generally term squeezed states those that have been transformed by a particular kind of 
unitary evolution named the squeezing transformation. 

To effect such an evolution in the laboratory requires the use of nonlinear optical elements and optical 
pumping 1571 l58l . The simplest example of an squeezing evolution involves the use of a nonlinear medium 
down-converting photons of a given frequency to two photons of half this frequency 1431 [581 . Let a be the 
annihilation operator for the mode to be squeezed, of frequency huj. Let b be the annihilation operator for 
an intense coherent field; the pump, of frequency 2hio. Let S'^^ be the strength of the coupling between the 
two modes, and the nonlinear coefficient of the medium inside which the squeezing evolution occurs. For a 
lossless setup the Hamiltonian for the interaction is given by 



Tri (pi,3 1 ) 1 ( X W 1 1 ) = P (x W )p3 W ) 



(1.33) 



P3 = 




(1.34) 



(1.35) 



For each photon of the pump annihilated, two photons of the mode are created, and viceversa. Given the high 

"^According to the Copenhagen Intepretation of Quantum Mechanics. 
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intensity of the pump field and its quasi-classical character, it's annihilation operator can be approximated by 
a complex amplitude, b k, (3. 

The unitary evolution operator, U{t) = e^^^ * for the time independent Hamiltonian in eq. ( 11.351 ) 
will therefore be 

5(0 = eHc*'i'-C'i'') (1.36) 

with ( = S'^' /? t, the phase of argument of the operator being that of the pump field. This evolution operator 
is unitary, therefore S{~C) '5^(0 = §"HC)- 

S{C) has been named the squeezing operator 1431 1471 . having a complex argument C = \r\ e**^. The 
modulus Ifl of this argument, |r|, usually called the squeezing coefficient is a characteristic interaction time 
and a logarithmic measure of the degree of squeezing. The squeezing phase if, equal to that of the pump field, 
will determine the phase-space orientation of the squeezing transformation. 

Using the ordering identities of eqs. ( 11.5b and ( 11.151 ). the results of the squeezing transformation on 
annihilation and creation operators, and on the displacement operator can be derived, 

§(C)a5t(C) =a cosh(r) - sinh(r) e*"^ 
SiOa^S^iC) cosh(r) - a sinh(r) e"*"^ 
S{OD{a)S\0 ^D{a cosh(r) - a* sinh(r) e'"^ ) (1.37) 

Given eqs. ( 11.311 ) and ( 11.371 ). the result of the squeezing transformation of the generalized quadrature i 
can be calculated; 

S[0 S\0 = -e" sin($ - ip/2) i('^/2+-/2) + ^-r ^^^^^ _ ^j^) ff'^/^) (1.38) 

Taking $ = (^/2; the expression in eq. ( 11.381 ) is simplified as the phase orientation of the quadrature x'-*-' 
coincides with the phase orientation of the squeezing operator. Thus we have, 

5(re2'*) x(*+-/2) s^re"'^) =e^x^^+-/^) 
A(5(re2''*) £(*+-/2) 5t(^e2.*))2 = g^"^ A(x(*+-/2))2 39) 

For r > 0, the variance of i*^*^ is scaled by a factor of e^^^ and the associated wavefunction is made 
narrower; this is termed squeezing in the language of quantum optics. While the opposite happens to the 
conjugate quadrature "expanded" by the inverse factor e^"^. Thus the squeezing transformation 

preserves the product of variances in Heisenberg's uncertainty relation, eq. ( 11.111 ). A quantum state for which 
that product is equal to the lower bound of j^, a minimum uncertainty state, will continue to be of minimum 
uncertainty after being "squeezed". 

With the purpose of simplifying calculations, we will take the squeezing operator's argument to be real; 
C, = r, where r can be negative. This is equivalent to the choice made for eq. ( |1.39t of squeezing phase- 
orientation; on quadratures a;(°) and p = x'''^/'^\ 

'We will usually refer to a signed, real quantity r instead of the non-negative modulus \r\. 
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The states usually associated with the squeezing operator S{C) are the coherent squeezed states Il47ll59l . 
produced by the squeezing transformation of a vacuum state; on which is then performed a displacement 
operation. Consider the case where <1> = = 0; 

|a;r) = S'(r)lO) (1.40) 

The coherent squeezed states are minimum uncertainty states like the coherent states, with the variance 
Ai^ = e~^''/4 reduced and A(p^ = e^'"/4 increased when r > 0. When r < 0, which corresponds to 
ip = TT, the p quadrature is the one "squeezed" and it's variance is reduced. The wave function in either 
basis of representation is Gaussian, and the averages of position and momentum are, respectively, the real 
and imaginary part of a. 

The infinite squeezing (and coupling x intensity of pump x time, as r = |/3|C'-'^-'t) limit for the squeezed 
states S{re^'^) are the x^^fl"^^ quadrature eigenstates. In such a way a physically feasible approximation 
for a quadrature eigenstate can be obtained by appropriate squeezing of an initial coherent state where a is 
wholly real (position eigenstate, r > 0) or imaginary (momentum eigenstate, r < 0), the limitation being 
technological. 

Quantum state preparation approximating the maximally entangled EPR states is possible given strong 
nonlinear interactions, control of the phase of the individual modes, and beam-splitters (see eqs. ( 11.251 ), ( 11.261 ) 
and ( 11.281 )). Two coherent squeezed states, squeezed in position and in momentum can be mixed by a beam- 
splitter obtaining a Gaussian state "squeezed" in nonlocal quadratures. Which, in the limit r co becomes 
an EPR state. 

The operation described above can be represented by the Two-mode squeezing operator BtI . Two-mode 
squeezing is the transformation preparing an entangled, symmetric Gaussian state state; from a two mode 
vacuum state vector | 0)_1, (8i | 0)2: 

§12(0 = e-'^'^^'^^+'^°i°= C,= re"^ (1.41) 

It is straightforward to show that the two-mode squeezing operator can be written as the product of two 
squeezing operators for different modes and a beam-splitter transformation; 

^12 (r) = Bi2(V4) 5i(r) S^i^-r) (1.42) 

where the two squeezing operations represent the nonlinear interaction (see eqs. ( 11.351 ) and ( |1.36l l) trans- 
forming each of the initial vacuums into states squeezed, respectively, in the quadratures x\ and p2. The 
symmetric beam-splitter transformation mixes the two squeezed states (see eq. ( 11.281 ) and discussion above) 
into an approximate EPR state, the two-mode squeezed vacuum; |C)as- 

Recalling the beam-splitter transformation's operation on a two-mode quantum state wavefunction, it 
is easily seen that the two mode-squeezed vacuum is Gaussian in wavefunction. The two-mode squeezed 
vacuum is an entangled state, and correlations arise in the measurement of observables other than the nonlocal 
quadratures. In the Fock state basis representation; 



oc 

|f)i2= (cosh(r))"^ ^ (tanh(r) )" |n)2 

n=0 



(1.43) 
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The state vector is symmetric under exchange of modes 1 and 2, having an even overall number of 
photons. A measurement of photon number must give the same result for modes 1 and 2. For the two-mode 
squeezed state ( rii — ^2) =0, even if the overall photon number rii + 712 is not unambiguously defined. 

The factor (tanh(7'))" in eq. ( 11.431 1 makes for a probability of measuring (overall) 2 n photons that 
decreases exponentially with 2n; and ensures the convergence of the infinite sum for finite r. The infinite 
squeezing limit r ^ cx) of eq. ( |1.43l l is the EPR state of eq. ( 11.281 ). In this limit the infinite sum does not 
converge and the normalization constant (coshr)^^ becomes infinitesimal. 

To study the effect of two-mode squeezing on non-Gaussian states' wave functions and phase-space 
representations; it suffices to study the transformation of modes' operators ai and 0.2', 

SuiO a^ SUiO = cosh(r) a, ~ e"^ sinh(r) a] for z ^ j = 1, 2 (1.44) 

a Bogoliubov transformation mixing mode operators. 

Analytic functions of the operators 0,1^2, o\ 2 will be transformed in an analogous manner to that of 
eq. ( II.25I 1; with the Bogoliubov transformation of eq. d 1 .441 ) substituting for the beam-splitter transformation. 
An interesting and useful (for the purposes of this work) example is that of a product of displacement 
operators in two modes, 

^12(0 D2{a2) SUO = i5iK) 52(4) (1-45) 

a- = cosh{r)a, + e*"^ sinh(r) a* for i ^ j = 1,2 (1-46) 

1.3 The Wigner and the Wigner characteristic function 

In this section, we introduce phase-space representations; of quantum states' density matrices, and of any 
square-integrable operator ll46ll . In particular, we introduce the Wigner function II6OII and the Wigner char- 
acteristic function 1461 as particularly useful for the representation of CV systems because they are explicit, 
analytic functions of phase-space "coordinates" that transform in the same manner as the wave functions of 
quantum states BSl . Though the Wigner function is not a genuine probability distribution in momentum 
and position, the two conjugate functions allow the calculation of quantum mechanical averages, including 
partial traces, to take the form of integrals over the complex plane of phase-space. Unitary evolutions and 
measurements over a multi-mode quantum state likewise take the form of simple unitary transformations on 
the arguments of the functions, and projections on adequate eigenstates. 

1.3.1 Phase-space representations: mainly Wigner function and characteristic func- 
tion 

The Wigner function was proposed and chosen to be an Hermitian form (real scalar) of the density operator 
of a quantum state fulfilling a number of desirable conditions; to be real and bounded, to transform according 
to the same rules for a classical distribution of probability; to produce the quantum mechanical averages 
pertaining to the density operator Bsl ; to give the appropriate probability distributions when integrated. 
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Such a function was proposed in ref. ||6()I in the form 

W{a = x + ip)= J dy (x-^\p\x+^) e^'Py (1.47) 

with the unit system defined so h = i. The Wigner function is real, and can be demonstrated to exist and 
be square-integrable for any density operator P6l . It has been termed a pseudo-distribution and as such, a 
conjugate function, it's Fourier transform, the Wigner characteristic function has been defined 



(1.48) 



W{a) = vr-i J d^i e"«*-"'« xiO (1-49) 

with ^ = w + iz being a conjugate phase-space "coordinate". The characteristic function for a density 
operator, being the Fourier transform of the analytic, real Wigner's function is thus an analytic function. 
We remark here that "Wigner" representations of suitable operators other than the density operators can be 
calculated, and are necessary to the calculation of physicalle relevant quantum mechanical averages using this 
representation. Let us review the forms (in a strict mathematical sense) corresponding to quantum mechanical 
averages and normalization conditions. 

The normaUzation condition on density operators; Tr(/5) = 1 is equivalent, in the Wigner and Character- 
istic Function descriptions to 

Tr(p) = TT-i J d^ i XpiO - j Wp{a) = 1 (1.50) 

Finite Wigner function and characteristic functions can be derived, and used to calculate quantum me- 
chanical averages for any operator fulfilhng the condition that the Hilbert-Schmidt norm Tr(FtF) be finite, m 
other words that the operator be bounded 1461 . The density operators are, furthermore, trace-class operators 
with trace equal to 1; therefore Tr(/5^) < 1, and 

Tr(p2) = j d^^\xp{0? = j d^a\Wp{a)\' < 1 (1.51) 

The quantity Tr(/5^) is a quantitative measure of the purity of the quantum state, and is thus named. 
Quantum mechanical averages or expectation values of the product of at least one trace-class operator 
and a bounded operator are finite and are given by the form Tr(FG'); 



Tr(FG) =^-1 J d'CxFiO XG(-e) (1.52) 
= 7r-i I d^ a Wpia) Wcia) (1.53) 



where Wp, Wg and xf, Xg are the Wigner and characteristic functions corresponding to F, G respectively. 

The Wigner function can be easily shown not to be a genuine probability distribution for density operators, 
even though it produces expectation values and is normalized to 1 Uke a classical probability distribution. Let 
/5$ = 1^) (^1 and = j^") (^Pl, density operators for two pure states. If the states are orthogonal, then 
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( $ I 4* ) = Tr(/5$ pq,) = 0. The integral of the product of two nonzero Wigner functions in eq. ( 11.531 ) is 
equal to zero. Hence, the Wigner function of at least one of the density operators must have negative values. 
The Wigner function cannot be a genuine probability distribution for all quantum states, and is thus termed a 
pseudo-probability. Conversely, quantum states with a positive Wigner function that is a genuine probabihty 
distribution cannot be orthogonal with each other. 

Quantum states of a nonclassical character have a Wigner function IfTTl [T3l with negative values in parts 
of its domain; the negative volume of the Wigner function has been proposed as a measure of the nonclassical 
character of a quantum state IMI . 

The Wigner characteristic function is the trace of the product of the density operator and the displacement 
operator ll46l 1481 : the mean value of the displacement operator for the density p. Recall that displacement 
operators form an orthogonal basis for the representation of square-integrable operators (see eq. ( I1.23l l and 
that any bounded operator may be represented by a linear combination of displacement operators P6l . The 
characteristic function is the coefficient for this representation, and is necessarily square-integrable (see 
eq. (I1.51I I) when the operator itself is trace-class. The characteristic function is also given by 

x(0=Tr(p5(0) (1.54) 

which for a pure state p = (V'l becomes 

X(0= (1.55) 



Generalizing the characteristic function to more than one mode is straightforward; the product of op- 
erators to be traced over in eq. ( 11.54b being that of the density operator and the commuting displacement 
operators -Di(^i) 02(^2) ■ ■ ■ Dn{£^n), for modes. 

A most important property of the characteristic function, is that transformations on density operators of 
quantum states such as those we have reviewed in the preceding section will correspond, in the language 
of the characteristic functions, to transformations on the arguments of the displacement operators (see 
eqs. ( 11.27b . dl.lSb . ( 11.37b and ( 11.46b ); which are the arguments of the characteristic function. 

The Wigner function and Wigner characteristic function are associated to a symmetric ordering of op- 
erators 14611481 ; where ordering means the left to right ordering of the non-commuting a and operators 
inside other operators such as density matrices and observables. For example: the number operator n = a^a 
is symmetrically ordered as 2^^{a^a + aa^), normally ordered as a^a, anti-normally ordered as aa^ . 

Other phase-space representations exist based on alternative orderings. In the manner discussed above, the 
-^(q;) lf38l representation and the Q{a) 1621 function are conjugate with characteristic functions associated 
with the normal and anti-normal ordering of operators, respectively. Particularly, with the ordering of the 
displacement operator of which the characteristic function is a coefficient of representation; 

{3(0} =e^''-^'' 
5(e,s)=e^l«'l/2{5(a)} 

5(^,-1) =e-«*"e«°' (1.56) 
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where the index takes the values s = 1 for normal ordering, s = for symmetric ordering, also denoted 
by the brackets {...}, and s = — 1 for antinormal ordering. It is obvious from inspection of eq. ( 11.541 ) and 
eq. (II. 56b that the transformation between characteristic functions associated with different orderings is just 
a power of e'^^'/^. 

The characteristic function of the Wigner distribution, similarly to a classical characteristic function, is 
the generating function for the symmetrically ordered moments of annihilation and creation operators BTI : 



simpUfying the calculation of statistical moments; specially quantum mechanical averages and covariances. 



(1.57) 

? = C*=o 



1.3.2 Characteristic functions of the two-mode squeezed vacuum and EPR states 

The characteristic function for a Gaussian, two-mode squeezed vacuum state |0)a (8) |0)h is given by 

X5(Ca,6) =e-V2(lc;i^ + l?^l^) (1.58) 

where the relation of the pair of variables , to the pair of variables , Cfc is the Bogoliubov transformation 
described in eq. ( 11.461 ). Other, displaced two-mode coherent states have a characteristic function similar to 
the one in eq. ( 11.581 ) with a multiplying phase factor. For an overall phase-space displacement of, say x' + ip' 
on A mode this factor is given by e^''^° ^ - p ) 

The infinite squeezing Umit of the Gaussian two-mode squeezed state is the EPR state. In a CV setting, it 
is the maximally entangled state 1451 of two nonlocal quadratures Xu,pv, as defined by the unitary transfor- 
mation of eq. ( 11.271 1. Given the wave function (see eq. (11.28b ). the Wigner function of the state can be derived 
easily from it's definition in eq. ( 11.47b . 

Wx'^,p'^{aa, ab) = C S{cos{Q)xa - siii(e)x6 - x'^) 6{sm{Q)pa + cos{Q)pb - p'^) (1.59) 

where C is a normalization factor that becomes infinitesimally small as r oo. The characteristic function 
of the EPR state is easy to calculate as the Fourier transform (eq. ( 11.48b ) of eq. (11. 59b ; 

X<.P'„(ea,6) - C S{za sin(e) + Zb cos(e)) S{wa 008(8) - Wb sin(e)) e2*(^°< -c(e)-™„p; csc(e)) 

(1.60) 

which is, again, not square-integrable. However, the product of this characteristic function and a square- 
integrable characteristic function can be traced over (integrated) in the manner of eq. ( 11.52b . 



1.4 Universal quantum teleportation (in Continuous Variables) 

This section will review both the universal protocol for Quantum Teleportation of a quantum state (and apply 
it to a CV setting) and the fidelity of teleportation in the language of Wigner characteristic functions. 
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1.4.1 Quantum teleportation 

The quantum teleportation protocol involves the transcription of the unknown quantum state of a physical 
system in (named input) onto the quantum state of another, similar physical system 50 that is remote with 
respect to in. The protocol itself is a projective measurement of the maximally entangled states for that 
setting Il25ll26ll63]| . the Bell states (in Continuous Variables the EPR states) on the joint system consisting 
on state in and another, similar system 

Quantum teleportation is not "cloning" the quantum state, as this operation is impossible for arbitrary 
states |l64| . It is no scheme of quantum measurements on the input state or of repeated measurements on 
identical instances of preparation as is done in quantum tomography 16511661 . The state itself is unknown and 
is to remain unknown through teleportation; in principle there is only one copy available for teleportation. 

The transcription of an unknown state from in to physically separate system B by measurements on 
systems A and in is only possible if a quantum correlation exists between the states of the two systems A 
and B (the joint, entangled state of both being named resource) and if a correlation is forced on the A, in 
state by interaction between the two systems before measurement; there is no Bell state measurement without 
previous interaction ll67l . 

An obvious maximally entangled or Bell state for a quantum system Hilbert space H of dimensionality 
d can be formulated with ease 1261 . Let {|n)} be a complete basis of representation of H; the simplest Bell 
state has a state vector 

1*0 = d-'/'' ^ e-*^" \ n)A \n)B (1.61) 

n 

belonging to the joint Hilbert space Ti. (E) H. This state has the same form when written in terms of any 
complete basis of representation of H. Let the basis { |n) } be made of eigenstates of an observable operator n; 
Alice 's measurement of will determine the result for Bob 's measurement of fiB- It follows, on inspection, 
that the Bell state of eq. (II.6II 1 is an eigenstate of nonlocal observable fiA — fiB with eigenvalue 0. nonlocal 
observables having as eigenstates the maximally entangled Bell states are termed Bell observables; their 
measurement is taken to be a projective measurement of the Bell state Note that the squeezed vacuum 
of eq. ( II.43I 1 and its infinite squeezing Umit, the EPR state of CV systems are eigenstates of a nonlocal 
observable (the difference in photon number of the two systems). 

Other maximally entangled eigenstates, associated to other eigenvalues of Bell observables can be pro- 
duced by a local unitary transformation on eq. (11.61b ; 

\^a)AB^ UA{g)\^o)AB (1.62) 

In this case Alice has performed the operation U{g). To produce all the possible Bell states, U[g) must 
belong to a group of unitary transformations {£/((?)} having an irreducible unitary representation in the space 
of operators acting on 7i. Schur's lemma for the density operators of the Bell states takes the form: 

j dg UA{g) I *o )ab ( *o \ab U\{g) = d"! Xa ® is (1.63) 

^belonging to a party named Bob. 

^belonging to a party named Alice, also in possession of the in system. 
^Following the Copenhagen Interpretation of Quantum Mechanics. 
'Also termed infinite dimensional 
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with the integration over the argument g being replaced by a sum whenever the group is discrete. The factor 
becomes infinitesimal for CV systems. 

The identity in eq. ( 11.63b is a proof that the group consisting of density operators of Bell states form a 
Positive Operator Valued Measure ll25l (abbreviated POVM). It is also a proof that the Bell states form a 
complete, orthonormal basis for the Hilbert space TC <E)'H; eq. ( |1.63l l is a completeness relationship for Bell 
basis states. 

For the 2-dimensional (spin 1/2) Hilbert space of the original teleportation proposal lISTl . the unitary 
transformation group consists in the set of the Pauli operators { 1 , }. The Bell basis is 

given by 

|T>a|T)b + e*U)AU>B 

I T)aU)b + e''^U)A| T>i3 (1.64) 



with = 0, TT. The Bell observables are given by 



(^) 



(1.65) 



In a CV setting a 2-dimensional Hilbert space (of bosonic excitations, though) can be defined by a 
truncation of the Fock basis to the first two states {|0), |1)}. The most general superposition state for this 
truncated space is of the form cos(£) | ) + e*® sin(£) 1 1 ). The Bell observables are fiA — nB and fiA + nB 
and their eigenstates have a form similar to that of eq. ( 11.64b . 

For a full CV setting the group of unitary transformations producing the Bell states is that of the dis- 
placement operators and the Bell basis is made of the EPR states of eq. ( |1.28b . The state in eq. ( |1.61b being 
the infinitely squeezed two-mode vacuum (see eq. (11.43b ): having nonlocal quadratures' (Bell observables) 
eigenvalues pa+ Pb = and xa — xb =0. The displacement operation, when performed by either Alice 
or Bob will produce all the other possible values for the nonlocal quadratures. The POVM for the Bell-basis 
measurement of Bell operators will be the homodyne measurement POVM ||63l . 

For the quantum teleportation protocol, we have a three-party joint system of the entangled resource AB, 
and the input in with a density operator pin ® pab- The teleportation protocol consists in the measurement of 
the Bell observables (projective measurement of the Bell states) over the modes A, in in Alice 's possession. 

Given the POVM defined in eq. ( |1.63b . a measurement with a definite result g would result in a state, after 
measurement (see eqs. (11.32b . (11.33b and (11.34b and the associated discussion of projective measurements); 



(1.66) 

which is generally not normalized. To simplify and restrict exposition to Bell resource teleportation, let 
the input state be the pure state \ip)in', let the resource be |'i'o)yiB of eq. (11.61b : lastly, let the basis of 
representation for the states { | n ) } be orthonormal. The outcome of the projective measurement of eq. ( 11.66b 
with a result g is the state 

d-i J2 e'*"'-**" {n'lAinUnULig) \n) A\n) b \ ^)^n = d'^ ^ ( n | [/t(g) | ) \n)B (1.67) 

n,n' n 
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Given the complete and orthonormal nature of the basis {| n )}; the state in eq. ( 11.671 ) is the original input 
state (now "living" in system B instead of in), on which the transformation (g) has been performed. 

Note that the Bell state for the outcome of measurement g and the resource state are related by the 
transformation \'i>g)AB = UA{g)\'^o)AB- Moreover, the generalization to a teleportation protocol where the 
resource state is a Bell state different from l^'o)^!^ is straightforward; as {U{g)} is a group, with a group 
composition law: 

U{g')U{9") = C/(g'..g")e^*^^''^"^ (1-68) 

where e'*^^ Ms a phase factor (f.e. the phase in eq. d 1.181 ) for displacements); and the operation indicated 
by g' ■ g" is algebraic. The phase factor and the arguments g under the operation • are associative, with an 
inverse element and an identity element. Thus, the transformations on a Bell state not equal to |^'o)yiB can 
be compounded with ease. 

To finish the teleportation protocol. Bob must apply the transformation UB{g) onto the state B in his 
possession. The result of measurement is random in ?i fundamental manner; to know which transformation to 
apply Bob has to receive a communication from AZ/ce via a classical communication channel. 

Owing to the random nature of projective measurement. Bob 's output state should in reality be an ensem- 
ble of all the possible outputs corresponding to possible results of such measurement (as in eq. ( 11.341 )). We 
have explained teleportation with a Bell state resource, entailing perfect transcription of the input state for 
every measurement result. For this special case, the ensemble nature or "mixedness" of the output induced 
by teleportation does not exist as all results produce the same output. In a CV setting, and in any other 
setting where the protocol does not rely on a Bell state resource and produces different outputs according to 
different measurement results, the mixed nature of the output state is relevant and is to be taken into account. 
In refs. 13311681 and 1691 (dealing with CV teleportation) integration over the possible outcomes of homodyne 
measurement is performed on the un-normalized output state of eq. ( 11.661 )). 

In the CV protocol l33l Alice communicates the results of homodyne measurements of nonlocal quadra- 
tures to Bob; he applies this same displacement to system B in his possession. For a review of the CV 
teleportation protocol as a projective measurement of pure state wave functions; see also refs. iTOl [7111721 . 

1.4.2 The fidelity of teleportation for Continuous Variables 

The most widely accepted definition l39l for the fidelity between two quantum states; pout, the output state 
of quantum teleportation and pi„, the input state is given by the trace of the product of both density operators. 
In the the phase-space representations used in this work, such a fidelity is given by 



=Tr{pin Pout) (1-69) 
= 7r-i y" d^a WM Woutia) (1.70) 

= n-' J d^^X^n{OXout{-0 (1-71) 

The positivity and normalization properties of all density operators (see eq. ( 11.51b ) determine that the 
fidelity is bounded from above and below: < J- < 1. Supposing that at least one of the states is pure; say. 
Pin (see ref. l39l ). a fidelity of 1 means that /5j„ = Pout- A fidelity of would mean that the two states are 
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orthogonal. 



Generally, the "practical" lower bound to fidelity is taken to be the classical teleportation threshold or 
classical fidelity, for a given input and a separable two-party resource; considered a worst, limit case of a 
class of possible entangled resources. For the classical teleportation of an arbitrary coherent state by separate 
two-mode vacuum states [39|, it is Tds = ^■ 

Such a classical threshold can be calculated for a given input, and a class of resource states. First, define 
a separable, "limit" resource state, with density operator pA ® Pb and an input state pm (or an ensemble of 
likely inputs). To perform the teleportation protocol a projective measurement of a Bell state is performed by 
Alice on rhoA ® Pin- Probabihty for result g given by 



The ensemble state of all the possible "corrections" performed by Bob on state in his possession will be 
the output state for the teleportation protocol; 



Pclassig) = Tr(/5A ® Pin \ )in,A ( *g |m,yl) 



(1.72) 




(1.73) 



a 



The classical fidelity is given by eq. (|1.69l i; calculated for the input state pin and the output state of 
eq. ( 11.73b . This fidelity can be modulated and optimized by changing constraints on the hkely input states' 
ensemble and on the character of the limit, separable resource. 
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Characteristic function formalism for 
CV teleportation 

The universal procedure for teleportation (see section 11.41 ) of input quantum state (designated in) using an 
entangled resource shared by parties Alice and Bob (designated A , B) consists on a projective measurement 
by Alice of A and in joint system onto a Bell state. Followed by a local unitary operation on B by Bob, based 
on measurement results communicated by Alice. 

In this chapter, we will describe quantum teleportation using characteristic functions as a means of 
representation. In this representation it is straightforward to obtain an elegant expression for the teleported 
state; for all (one-mode) input states and all (two-mode) resource states in a CV context, and to introduce 
simple modifications of the teleportation protocol which obtain general, elegant expressions for the output 
state. For these reasons, we speak of a characteristic function formalism for CV teleportation. 

Mixing in external modes by means of beam-sphtters, making projective measurements on the modes 
arising from such beam-splitters, slightly modifying the nonlocal measurement that constitutes the basis of 
quantum teleportation are a few examples of the possible modifications on the basic CV teleportation protocol 
that can be studied with ease with the characteristic function formalism. 

2.1 Teleportation in the language of Wigner functions 

The maximally entangled states in a CV setting are the EPR states see (eq. ( 11.28b ). These are the eigenstates 
of the nonlocal quadratures Xu andpv defined in eq. ( 11.27b . Entangled states of these quadratures are achieved 
in CV settings from separable states by means of beam-splitters (as illustrated in section [r.2.2l ). 

To perform quantum teleportation in this setting 133)1 . Alice must produce an entangled state of modes A 
and in by means of a beam-splitter; then she is able to perform measurements of nonlocal quadratures ll67l 
by homodyne detection. The experimental setup in fig. ( 12.1b illustrates such a scheme. 

Initially, the joint system of input and resource is in a state given by the Wigner function 

WAsiaA-, aB)W,n{a,n) (2.1) 

' Mostly related to realistic, noisy measurement or preparation of the resource 
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Entangled Resource 
Generation 



Figure 2.1: The Experimental Setup for Continuous Variable Teleportation. Alice mixes input mode in with 
mode A of the entangled state that is used as a teleportation resource by means of the beam-splitter "0". 
An homodyne measurement of nonlocal, commuting Xu , Pv quadratures is performed, "collapsing" the joint 
state of modes A and in into an EPR state and transferring the initial state of mode in onto mode B; save for a 
phase-space displacement corresponding to the actual values of the quadratures measured. This displacement 
is corrected for by Bob on receiving a communication of the measurement results from Alice. 



The "(?" beam-splitter, with transmission coefficient cos(6') (fig. (12.1b ) mixes the modes incoming modes 
A and in, producing the outcoming modes u and v corresponding to nonlocal phase-space variables a„ and 
a„. The state of the system after the beam-splitter operation is given by 

WABiaAiau,av), as) W:i„(ai„(a„, a^,)) (2.2) 

where variables a^(a„, a^) , ai„(a„, a,,) are given by the inverse of the transformation in eq. (I1.27l i (see 
also eq. (11.26b ). 

On such a entangled state of A and in modes homodyne measurements of the nonlocal, commuting 
quadratures a;„ and p„ are performed; thus a Bell observable measurement is realized. The state of part B 
of the system after a projective measurement with results x'^ and p'^ is represented by the partial trace over 
A and in (see eq. (11.53b of the Wigner function of the system (eq. (12.2b and the Wigner function Wx'^y^ 
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(eq. ( 11.591 )) of the respective element|§ of the Bell measurement; 

{ctB ; x'^,p'v) = TT ^ J (fain(faA WAB{aA{au,ay), as) Win {a in {a u, ay)) 

X C d{cos{6)x.in - sm{0)xA - x'^) S{sm{9)pin + cos{0)pA - p'y) (2.3) 

x'^,Py) in eq. (|23] i is not a normalized Wigner function (as is the case in eqs. ( |1.33b 
and (II.66I 1). It is the product of the probability P{x'^,p'^) of a measurement result x'^ , p^ and the Wigner 
function of the state of the B part of the system after a measurement giving such a result; given by the 
conditional pseudo-probability Wmeasured{o:B \ x'i^,p'v)- The probability can be obtained from eq. ( 12.3b by 
tracing out the B mode; 



(2.4) 



Perform the integral over ua in eq. (I2.3l l. yielding a convolution integral of the resource and the input 
states; 

7r"2C f 

WmeasuredictB] x'^,p'v) = , . , 77^ / Olin Win{ain) 

\sm{e)\ \cos(0)\ J 

xWAB{xinCot{9) - x'y^csc{e) + 1 {py sec{9) - pinta.n{d)) ; ob) (2.5) 

Up to this point, we have not discussed the operations that Bob must perform on B mode. Only to see 
clearly what these might be, we will consider an ideal resource, the EPR state of Wigner function 

WABiaA] ub) = 2C5{xA - xb) 5{pa + pb) (2.6) 

Substituting this resource in eq. (I2.5l and integrating yields the resulting, un-normaUzed Wigner function; 

2 7r"2^2 

-— — 7^W„ (xs tan(6') + < sec{9) + i {pb coi{0) + p',v) csc{e))) (2.7) 

I sm^t/jl I cos(6')| ^ 

where the probability of obtaining the arbitrary result x'^ , p'„ of measurement is easily seen to be the vanishing 
constant 2 C^TT-^ (|sin(6')| \ cos{e)\)-\ 

To recover the input state at his end. Bob must perform two unitary operations on the B mode; one of 
them after knowing the results of the measurement, communicated by Alice via a classical channel. Namely, 

• A displacement in phase-space of (a;^ sec(6') + ip[,csc{d)). For homodyne measurement apparatus 
with gain coefficients < < 1 and < < 1 for x„ and p„, respectively; this displacement 
will be {gxx'y^ sec{9) + i gpp'^ csc{6)). Thus the operation is D{{gxx'^ sec{9) + i gpp'y csc{9))). 

• A transformation xb + ipB — ^ o-xb + ia,~^PB, with a = coi{9). Therefore, is is a squeezing 
operation 5(r), where e^'' = cot{9). Obviously, cot{9) is part of the description of the beam-splitter 
inside our experimental setup and can be known beforehand by Alice and Bob. 



^Namely, the density operator of the EPR state with eigenvalues x'^ , p'^ . 
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The complete unitary operation applied by Bob would be given by the operator 

C/i3 =^B(-ln(cot(6l))) 5b(«scc(6I) + ip'^csc{e))) 

= 5b((4csc(0) + ip[,scc{e))) SBi-Hcom)) (2.8) 

The squeezing operation performed by Bob can be made unnecessary. Consider a teleportation resource 
of the form 

WAsiaA; as) = C5{xAt&TL{e) ~ XB)S{pACot{6) + pb) (2.9) 
instead of that of eq. ( 12.6b . This would give the end result 

I . ,nM I 7^^™(^B + <sec(0) +i{pB+ p[v)csc{e))) (2.10) 

thus, in this case. Bob must perform only the displacement above described to recover the input state and 
realize teleportation. 

Note that a squeezing transformation XA + ipA — ^ cot{9)xA + i tsin.{9)pA (identical to Bob 's squeezing 
"correction") on mode A of the resource in eq. ( |2.6l l would result in the resource of eq. (12. 9t . A displacement 
of {x'^ csc{9)+i p'^ sec{9)) on mode A of the resource in eq. ( 12. 9t would yield the EPR state Wx'^y^ , identical 
to the state of the A and in modes after the Bell measurement with results x'^ , p'„ . The unitary transformations 
on mode A described above; transforming the EPR state in eq. (I2.6l l into the EPR state Wx'^y^ of eq. ( I1.59l l 
are identical to those applied by Bob to mode B to realize teleportation (in eq. ( 12. 8t . 

We have just shown that teleportation using an asymmetric beam-splitter of angle 6 together with the 
simplest maximally entangled resource (eq. ( 12.6b . will result in added squeezing on the output state. That 
squeezing will have to be corrected by the application of the inverse squeezing transformation. Previously 
applying this inverse transformation to mode A of the resource will eliminate the need for such a correction 
by Bob. 

The use of an asymmetric beam-splitter for mixing the A and in modes would be desirable only if the 
teleportation protocol were intended to produce an squeezed output. The effect of having an asymmetric 
beam-splitter in the experimental setup can be mimicked by appropriate local transformations on the resource 
state. 

In order to simplify exposition, and assuming that we have no further use for additional squeezing of the 
output state we will take the beam-splitter "9" to be symmetric (6 = 7r/4) and without a phase. The use of a 
symmetric beam-splitter with a phase would only rotate in phase-space the displacement to be performed by 
Bob. 

Using a symmetric beam-splitter, and after Bob's correction (the displacement in eq. (12.8b ) we have for 
the output state of the system (see eq. (12.5b ): 



Wout{aB\ x'^,p'v) = 27r J Win{ain) 
X WABix^n - '2}l^x'^ + 1(2^2 _ p,„); - {'^^^'^Q^x'^ + iT}'^ g^pJ,)) (2.11) 

This output of teleportation is un-normalized, and so far dependent on the outcome x'^,p'^ of the Bell 
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measurement. It is the product of a probability for such an outcome and a conditional Wigner function (see 
eq. (I2.3l i. Given 

• the fundamental randomness of the results of quantum measurement, particularly Bell measurement of 



• that teleportation is, in principle performed in an absence of any knowledge (by Alice and Bob) of the 
input state; thus in the absence of any knowledge, even statistic, of measurement results. 

• that teleportation is performed in an automatic manner by Alice and Bob, without change to the 
experimental setup due to the knowledge of particular set of measurement results. 

• that for a realistic fidelity of teleportation, it is necessary to consider all the random outcomes of 
measurement, modulated by their probability; a fidelity coefficient depending on a single random result 
is not acceptable because it cannot be repeated reliably. 

• that a conditional output state, on a random result is not an acceptable answer for a teleportation output, 
as it comes about randomly. 

it becomes evident that a final output state that is acceptable is given by a mixture of conditional states (see 
eq. i2.3i ) with the probability of measurement in eq. (12. 4t . Therefore, integration over x'^^ and p'^ of eq. ( 12.1 lb 
will yield a normalized Wigner function that is an ensemble of conditional output states corresponding to 
individual measurement results. 



Wont [a b) ^ J dx'^dp[, WoutiaB; x'^,p'v) 

= 2n-^ J d^a.,n Wi„(a,„) 
X / dx'^dp', WABix^n " 2^^^ x'^ + I {2^/^ p', " K„); as - (g.2i/2< + igp2'/^p',)) (2.12) 



The Wigner function Wout{ctB) is the outcome of the convolution of the Wigner function of the input 
Win and a bipartite (on in and B modes) Wigner function given by 

A'(a.„; as) = 2 tt-^ ^ dx'^ dp', WAB{x^n - 2}'^ < + i {2^1^ p', - p.,n); aB ~ (ff. 2^'^ x'^ + ig^ 2^'^ p'^)) 

(2.13) 

named the teleportation kernel in the literature ll69ll73l . It is easily seen that the kernel is the Wigner function 
of an ensemble (with constant, flat probability) of Transfer Operators 12611741 (one for each value of x'^ , p^), 
having the Wigner function 

WAB{x,r. - 2^/2 + I (2^/2 - as - {g.2'/''x'^ + igp2'/^p'J) (2.14) 

Let the characteristic function XoMi(^s) be the Fourier transform of VFotit(as) (see eq. ( 11.481 1; 

Xouti^B) = 2 7r-3 y d^ae J d^a.^ J dx'^dp', ^M^^zs-pb^b) 
xT^™(a„) WAB{x^n - 2^/2 + i(2i/2p^ _ aB - {g.2'/^x'^ + igp2^^^p'^)) (2.15) 
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Making the substitution = = and multiplying the integrand in eq. ( 12.151 ) by 

g± 2 1 (g, ZB xZ - sp -^B P^') obtain 

X WAB{xi„ - a;" + - Pin); as - {gxx'l + igpp'^)) (2.16) 

This is an explicit Fourier transformation over the a^, x" + ip" variables; as well as a convolution over 
ain- Therefore, the characteristic function of the output state is straightforward to calculate: 

Xout{£.B) = XAsigpWB - igxZB; ^b) Xin{gpWB + igxzs) (2.17) 

A result equivalent to that obtained in ref. ifTSl for a symmetric beam-splitter using the transfer operator 
formalism. 

Keeping the beam-splitter asymmetric, and having Bob perform the squeezing operation described in 
eq. (12.8b will result in output state 

Xout{£.B) = XAB{gpWB-igxZB;t'A'a.{6)wB+i coi(0)zB) Xin{gp ia.n(e)wB+igx coi(0)zB) (2.18) 

2.2 Teleportation in the language of characteristic functions 

The derivation of the CV teleportation output made in the previous section will now be repeated in the 
characteristic function representation and will be shown to be much simpler. The experimental setup for 
teleportation is illustrated, as before, in fig. (12. 11 1. Let the initial state of the joint system A, B, in be described 
by their characteristic functions 

Xab{^a; ^b) Xj«(6«) (2-19) 

There is a (symmetric) beam-splitter transformation from this initial state into the modes u and v and their 
associated quadratures, which, it can be seen easily (see eq. ( 11.48b ). transforms the conjugate phase-space 
variables £,in,£,A in a similar manner to eq. (11.27b to xi^. We will consider the EPR state's characteristic 
function (eq. (11.60b ) in terms of the variables £,in,^A- And realize the partial trace, or projection into the 
POVM element (for results x'^,p'^) thus, 

Xmeasured{£.B; x'^.p'v) = TT^^ j (f^incf^A XABi^A] ^b) Xini^in) Xx'^,p'^{~Cin , -^a) 
=7r"^ J d^iind^^A XAb{£.A] £.b) Xm(Cm) 

X 2C<5K„ - WA)5{z^n + z^)e-2'(--2V=<-"'-2V^pO (2.20) 

This characteristic function is, like its Wigner function equivalent in eq. (12.3b . un-normalized. It is the 
product of the probability of measurement 'P{x'.^,p'^) of a result x^,p'„ and the conditional characteristic 
function of the system Xmeasuredi^B\x'u, Pv) , on the aforementioned results. Tracing out the B mode in 
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eq. (12.21b will give us the probability of measurement; 

Pi^'u^Pv) = TT"^ J Xnieasuredi^B] x'^ , p'v) (2-21) 

With the purpose of having a look into Bob 's part in the teleportation protocol, we will consider the 
resource to be in a simple EPR state of the form 

Xab{£.a; £.b) = 2C S(wA - wb) S{za + zb) (2.22) 

Performing integration of eq. ( 12.201 1 with the aforementioned resource yields 

Xmeasurediis-, p^) = 4 Tf-^ X™ (^b) 6 ^ - <) (2.23) 



which is a product of the characteristic function of the input state (in mode B), with an additional 
phase-space displacement of — — i2^/'^p'^\ and a constant probability for every single result of mea- 
surement of A-K^^C"^. This result is consistent with that of eq. (I2.7l i (for = 7r/4). Thus, Bob must perform 
that which, to his knowledge (given the non-unit gain of the apparatus) is the opposite displacement operation, 
Dsigx 2^/^ x'^ + igp "^^^^p'v) to recover the input state; 

Xouti^B-, = e-2'Ksp2V^K-.,g,2V^<) ^'^^p'J (2.24) 

For the same reasons and on the same considerations exposed in the previous section; an ouput state 
conditional on a random measurement result is not acceptable, while an ensemble of conditional states with 
adequate probabihty is an acceptable output state. Therefore, integrate eq. ( 12.24b over x'^ and to obtain 
the normalized, ensemble state 

Xout{£.B) =C^^ j dx'^dp'^ Xouti^B; x'^,p'J 

= 2tt-^ Jdx'^ dp[, J d^^^n d^U e- 2 ' (""^ 2'''^'" - »^ 

X XAb(?A; £,b) Xin{S,in) (2.25) 

The integration of eq. ( 12.25b is entirely straightforward, giving the output state 

Xout{£.B) = XinigpWB + igxZB) XABigpWB - igxZB] £.b) (2.26) 

which is identical to that obtained in eq. (12.17b . after a much shorter and more elegant calculation. 

The formalism just outlined is general for any combination of resource and input states and gives a very 
simple expression for the output of Quantum Teleportation in CV. It is possible to use resources that are mixed 
states, reflecting the results of a conditional operation performed during the preparation of said resource. Or 
to construct input states that are mixtures of the states (or likely superpositions thereof) used to encode qudits 
in quantum information processing with adequate probabilities; thus constructing the general input ensemble 
of a quantum CV channel of teleportation, for which the fidelity can be calculated. 
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A basic analysis of eq. ( 12.261 1 shows that obtaining a result other than a random, "classical" guess 
depends on having a high measurement gain (g^ ,gp^ 1) and on having a resource characteristic function 
XAsigpWB ~ igxZB't^B) that is nearly constant. This last condition is fulfilled by states approximating 
the EPR state of eq. (|2.22t . For example two-mode squeezed vacuums at high squeezing r; or other states 
showing great similarity with a two-mode squeezed vacuum at high squeezing. 

2.3 Teleportation: a projective measurement onto a mixed state 

We have produced a simple formalism for the representation of the elementary transformations and projective 
measurements performed in CV teleportation and produced a compact expression for the characteristic 
function of a general output state, for all resource states. It is conceivable that the first and most obvious 
change in the teleportation protocol involves the EPR state onto which we project to represent an homodyne 
measurement (see section[L4]and eq. ( 12.201 )). 

The first choice if we are interested in the introduction of "imperfect" homodyne measurements would 
be a suitable mixture of EPR states. What would a mixture imply? That we are still doing an homodyne 
measurement of the variables Xu , pv and projecting onto an EPR state. We do not know which state precisely, 
even if the apparatus for homodyne detection returns a result x'^J^,p[,. If the apparatus is imprecise (not 
damaged or lacking calibration) the distribution of outcomes will be centered on the values returned. We can 
define such a projecting state as the mixture 



where Xx_,p+ {£,in,^A) is the EPR state of eq. ( 11.601 ) (O = 7r/4) for eigenvalues x- andp+. The probability 
distribution P(x_,p+; x^,p',) is required to fulfill the following conditions if the state in eq. ( 12.28b is to 
represent an imprecise measuring apparatus; 




(2.27) 



(2.28) 




P(x_,p+; x'^,p[,) = P(a:_ - < , p+ - p'J 



X- = X, 



■u 



P+ = P: 



(2.29) 



namely to be normalized and "centered" around x'^ , p[ 



The mixture of EPR states of eq. ( 12.281 ) for such a probability distribution is given by 




X S{zin + ZA)S{Win - Wa) 



(2.30) 



Let us define the Fourier transform ofP{x,p), the characteristic function 




(2.31) 
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Using the definition of eq. ( 12.311 ). we can write the characteristic function in eq. ( 12.301 ) in a more elegant 
manner; 

Xra^M^n,U) = V(2^l^w,^ , 2^'^ z,^) 2C ^{z,^ + za)^{w,^ - w a) ^ ^v^.;, - 2 V^p^ ) ^2.32) 

To study the corrections to be made by Boh we substitute the mixture state eq. ( 12.321 ) into eq. ( I2.20l i in 
place of the EPR state Xx'^,p'^ \ and take the resource to be an EPR state (eq. (12.22b ). yielding an output state 

(2.33) 

We have the product of the constant probability for a given result equal to that in eq. ( 12.231 ): of the 
characteristic function P; and of the characteristic function of the input state, displaced in phase-space by 
— (a;(^2^/^ + ip^2^/^). Bob will try to correct for this displacement by applying at his end the displacement 
Dsigx 2^/^ x'^ + igp 2^^^p5j) on B mode. Again, the apparatus is assumed to have non-unit gains gx and gp. 

Proceeding in the same manner of section l2?2l performing the displacement just described on the charac- 
teristic function in eq. (I2.24l i. and taking the final output state to be an ensemble of outcomes corresponding 
to measurement results x'^ , p'„ (see eq. (I2.25l l) will result in the final output state 



(Cs) = P(-2^/^gpWs, -2^^'^gxZB) XinigpWs + igxZs) XAsigpWB - igxzs; Cs) (2.34) 



This is the same characteristic function of the output state in eq. (12.26b multiplied by the characteristic 
function P(-2V2gp , -2V2g^ zb). 

The output state of the previous section is thus "smeared" in phase-space. This is easily seen in the 
Wigner functions' language; as this product of characteristic functions is the Fourier transform of the con- 
volution integral between the probability P and the pseudo-probability Wi„ - For example, a nearly constant 
characteristic function P (the Fourier transform of a sharply peaked function P{x,p) ^ S{x) S{p)) will yield 
an output state in eq. (12.34b approximating the ideal output state of eq. ( 12.26b . 

The natural choice for a probability distribution P{x,p) intended to represent a measurement apparatus 
producing outcomes (as far as regards the projection caused by the measurement) with a random deviation 
from the "real", mean values x'^,p'^ is the Gaussian distribution; 

P(a;,p) =e''+''e-'="^^'-'='°p' (2.35) 
Piw, z) = e-''-""'^' - (2.36) 



producing a multiplying factor in eq. (12. 34b of 

/o-2r2 2\ /n-Ss 2 2\ 

e-(2e 2:f5jg~(2e g^w^) ^237) 



Reducing the measurement gains g^, gp might improve the output state, for input characteristic functions 
XinigpWB + igxZB) that become constant at a much slower rate than a Gaussian distribution when the 
variables gxZB , gpWB grow small, or for input characteristic functions with a wide support, the Fourier 
transforms of narrow Wigner functions. 
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2.4 Mixed teleportation resources: "superimposed" over noisy envi- 
ronments 

We have calculated the characteristic function for a mixture of EPR states (eq. ( |2.32t in section [273] In doing 
so, we have stated some properties of the mixture probability (eq. ( |2.29l l) that ensure consistency with our 
hypothesis of an imperfect measurement and greatly simplify the form of the output of teleportation when 
the projective measurement is done onto this mixture (eq. ( 12.341 ). 

In this section; we will derive the characteristic function of a teleportation resource that is a mixture of 
pure states in phase-space. Our first purpose is to obtain the teleportation output when such a mixture is 
used as a resource. Our second purpose is to obtain within our formahsm a well-known result for normally 
ordered phase-space pseudo-distributions. Namely the result of superimposing a state of the radiation field 
over another, preexisting state having a P{a) function 1351 [3811761 . It is generally understood that a given 
(pure) resource state will propagate in a spatial volume or be prepared ideally within the environment given by 
the pure vacuum |0 , ) state. In ref. 1381 it is spoken of the switching on/off of two sources in order; which 
is equivalent to preparing one state over another (or letting one state propagate over the same spatial volume 
of the) state we have deemed to be an initial environment. The superimposition we will consider in this 
section consists in preparing (or letting propagate) a two-mode teleportation resource state when the initial 
environment state is different from the vacuum and has a P{a) function. We will use environment states that 
are generally considered "noise" for our purposes; for example, separable, two-mode thermal states. 

Given the pure state pab^ let us define the phase-space mixture of states in an similar manner to that of 
eq. (|223; 

^(Mzx) ^ ^-2 j ^2^^ ^2^^ ^AiXa , Pa) P Bi^b , Pb) PAB.aaiat (2.38) 

where P^i and Pb are probability distributions that fulfill the conditions set forth in eq. ( |2.29t . The density 
operator 

PAB,a,,:,ai = Dsi'^b) D^iaa) PAB DA{aa) Dsiab) (2.39) 

is that of an ideal, pure state resource; displaced in phase-space by aa and ab. Thus, the density operator of 
eq. ( 12.381 ) is equivalent to that for two superimposed modes in ref. l38l . Note that the density operator for 
any phase-space mixture; for example the mixed EPR state (eq. (12.27b is of a similar form, as phase-space 
displacements transform one EPR state into another EPR state. 
The characteristic function of eq. ( 12.381 ) is thus given by 

Xab'^\£.A; £.b) =7r"^ J (fua (fab PA{Xa , Pa) P Bi^b , Pb) 

^ ^2^{zAXa-WAPa)^2i{zBX,-WBP,) ^^^i^^^.^ ) (2.40) 

where displacement operations correspond to phase factors in the characteristic function. Define the char- 
acteristic functions and Pb as the Fourier transforms of and Pb (see eq. (12.311 )). The characteristic 
function in eq. ( I2.40l i will have the factorized form 

Xab^'H^a; Cs) = PAiwA , za) PBiwB , zb) Cb) (2.41) 



^We use the term superimposition to avoid confusion with the fundamental Quantum Mechanics postulate and concept of 
superposition. 
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The characteristic function of eq. ( 12.411 ) can be made equivalent to the superimposition of the resource 
state pab over initial states in modes A and B. If the Wigner functions for such initial states can be substituted 
for the probability distributions Va and Ps in the mixture. This can be done only for Wigner functions 
fulfilling certain conditions; being finite, nonnegative and normalized. In other words, for Wigner functions 
that are genuine probability distributions. In refs. Il38l[35ll76l . the equivalent condition is set forth that the 
initial states possess a P{a) function representation. 

The Wigner function that is the Fourier transform of eq. ( 12.411 ) is given by 



^ab*^''(Q^a; ocb) ^ tt ^ (fuacfab WAsiaA - as - at) PAixa,Pa);'PBixb,Pb) (2.42) 



which is a convolution integral of the Wigner functions of the states. This relation is equivalent to that of the 
P{a) functions for "superposed" excitations in ref. 138)1 . 

The general mixed states that we have used as substitutes for the EPR states in projective measurement 
in section 12.31 can be thought of as superimpositions of a pure state on a "noisy" initial environment. The 
probabiUty distribution in eq. ( 12.271 ) is nonlocal; therefore it must be the nonnegative Wigner function of a 
nonlocal state. General, Gaussian entangled states, having nonnegative Wigner functions would be the first 
choice. This kind of superimposition can be interesting; but it falls outside of the scope of our work. We 
consider entangled states "pure" resources to be prepared; not noisy environments to "prepare over". 

We have chosen separate probabilities and P b instead of a joint probability in eq. ( 12.381 it is our 
purpose to represent preparation of a resource over "noisy" separable states and study the changes wrought 
by this new preparation on separability and teleportation fidelity, with respect to a pure state prepared over a 
two-mode vacuum. 

Thermal states on modes A and B constitute a noisy initial state over which the pure resource state can 
be prepared, resulting in a mixed resource state. Thermal states 1381 fulfill the requisite of having a positive 
Wigner function and having a P{a) representation. They are prime candidates for the representation of the 
simplest noisy environments. A single-mode thermal field of mean photon-number nxh has a density matrix 
in the Fock basis 



and a Gaussian l48l [38l Wigner function that is a genuine probability distribution. Which, along with its 
respective characteristic function, reads 



The characteristic functions of quantum states chosen to represent noise will appear as multiplicative 
factors in the output of teleportation. Substituting the resource of eq. ( I2.41l i in eq. ( 12.26b yields the outcome 

XoiA^b) = Xin{9pWB + igxZB) Xab{9pWb - igxZB] (,b)Pa{9pWb , -QxZb) Pb{wb , zb) (2.46) 




(2.43) 




(2.44) 



(2.45) 



This output has in common with eq. ( 12.34b the presence of a multiplicative factor; though of a different 
character. The mixture in that case is that of the "homodyne" projective measurement as "induced" by 
an imprecise apparatus. In this case the factor corresponds to a phase-space mixture that is the result of 
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preparation of the resource in a noisy environment. 



2.5 Losses and noise admission in the projective measurement 

In the previous sections we have used the characteristic functions language to produce an expression for 
teleportation outcome. We have introduced modifications to the protocol that imply, in our view, imprecisions 
in the homodyne measurement. We have formulated a "noisy", mixed resource state. 

A modification to the teleportation procedure is the next step. We will modify the teleportation protocol 
by the introduction of fictitious elements, with the initial purpose of modelling a realistic, lossy and noisy 
homodyne measurement. Setups introducing fictitious elements to homodyne detection were originally 
developed and adapted to the study of non ideal, "realistic" homodyne detection of single mode quadratures 
in the Q{a) function representation 1341 : and later independently developed for non ideal teleportation in a 
coherent state representation ifTTl . 

The realistic homodyne measurement will be realized by the addition of two beam-splitters of rather high 
transmittances cos^(^) and cos^(</j), where </> « « 1; and two "external" modes that will be mixed by 
these beam-splitters with the nonlocal modes (those originally to be measure: u and v in fig. ( 12. Il l) before the 
actual process of homodyne measurement takes place. These beam-splitters will also deduct a small loss of 
intensity to the original nonlocal modes. 

The setup for this lossy, noisy homodyne measurement is illustrated in fig. ( I2.5l l (see the left half of 
fig. ( 12. It for a comparison). The beam-splitters "^" and "(p" and modes dp and eo are the new elements in the 
homodyne setup. For the external modes a choice can be made of thermal states, vacuum states, or general 
one-mode Gaussian states including the aforementioned as special cases. 

The teleportation procedure will (as always) consist in the performance of a projective measurement of the 
quadratures Xu,Pv These are derived from the initial "0" quadratures by the beam-splitter transformations 
for the modes illustrated in fig. ( I2.5l l, 

/a„„\ /2-V2 -2-1/2 
[a, J 2-V2 2-V2 



cos((p) — sin(iy9) 
sin((/3) cos{(p) 








\ad) 




fay] 










cos((/)) — sin(^) 
sin(0) cos{(/)) 




(2.47) 



With these transformations at hand, it is easy to write the Wigner function of the appropriate EPR state 
for homodyne measurement (see eq. ( |1.59t ) having eigenvalues xj^ andp^; 

Wx'^,p[Xauo,ado] a^ojOfeo) =C S {cos{ip) - sin{ip) Xdo - x'J 

X S {cos{(f)) - sin{(f))peo - p[,) (2.48) 
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Homodyne 




Figure 2.2: The Experimental Setup for Realistic Homodyne Measurement in CV Teleportation. The external 
modes do and eq are mixed with the nonlocal modes uq and vq by the "ip" and "0" beam splitters. The 
resulting u and v modes undergo homodyne measurement and "projection" onto an EPR state. 

The characteristic function for the EPR state is obtained in a straightforward manner by Fourier transfor- 
mation of the above expression, and is given by 

X<,pJ,(CHo,Cdo; S.vo,£.eo) ^TT^C S{Wdo) S{z,J (5(w„J(5(ZeJ (5 (sin(i^) z„„ + cos{ip)zdo) 

X 6 (sin(0) wy„ + cos(</)) Weo ) e~ 2 » (^.o < - -^o p'. -^(-W) (2.49) 

The state of the joint system AB, in to be measured after the first, symmetric (7r/4) beam-splitter reads 

X^ni2-'/^U + 2-1/2^) XABi2-'^''U - '2'^'^U\ £.b) Xe.t. uiU) Xe.t. v{U) (2-50) 

where Xext. u and Xext. v are the characteristic functions of the two external modes do and eo, respectively. 

We follow the teleportation protocol as described in section l2?2l for the physical system in eq. (I2.50l i: using 
for the projective measurement the EPR state of eq. (12. 49b . As before, we have Bob perform a displacement 
of _Db(x^2^/^ + ip^2^/^) to correct the output state; only with the gains = .gp = 1; for the procedure we 
are performing already includes losses in intensity and is intended to model a lossy measurement. The output 
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state will be, as before, an ensemble of all possible measurement outcomes. Thus, 

X X/Ib(2""^/^Cuo ^ 2^^/^^i,o; £,b) Xext.uiUo) Xext.vi^eo) 

^ Xa;^ , (^Mo I Cdo i ^vq j S,eo) (2.51) 

After some (entirely straightforward) integration of the above expression, we obtain the output state 

Xout,n{^B) ^Xin{cos{<p) WB + 1 cos{ip) zb) XABicos{(f)) Wb - 1 coa{(p) Zb] C-b) 

{i2 ' sin((y9) Zb) Xext. {2^l'^sm{(j})wB) (2.52) 

The first obvious trait of this state is the multiplication by the characteristic functions of the external modes. 
These have arguments that are scaled by the (small) factors siii((p) and sin(^), while the characteristic 
functions of input and resource are scaled by factors (close to 1) cos((^) and cos((/)) which are the transmission 
coefficients for the beam-splitters and "1^9", respectively. For small angles and ip, eq. ( |2.52t approxi- 
mates eq. ( 12.261 ) for gains given as g^. = cos{lp) and gp = cos{(f>). In the limit case </> = 95 = eq. ( 12.521 ) 
equals eq. ( 12.261 ), as the characteristic functions of the external modes with argument are constant. 

In sections |23] and l24l we have multiplying characteristic functions of probability distributions (which 
can be also Wigner functions) as multiplying factors in teleportation output states. But there are two main 
differences. First, it is not required by the formalism we have used that Xext. u and Xext. v be the Fourier 
transforms of nonnegative Wigner functions. Second and most important, the multiplying factors in the 
arguments of said characteristic functions are siii(iy9) and sin(0): not the gains; but rather the complementary 
quantities of those which may be considered measurement gains (cos{ip) and cos(0)). Only in the case where 
the Wigner functions of the external modes are positive and (f) ^ ip = n / 4 (or g,j. = gp = 2^^/^) are the 
two cases entirely equivalent. These two differences illustrate the "lossy" nature of the teleportation protocol 
outlined here, as compared with the merely "noisy" protocols. 

There is an obvious choice of external modes' states to represent losses and added noise in a teleportation 
setup 1341 1771 . It is the general one-mode Gaussian vacuum states, which may be squeezed and of non- 
minimum uncertainty (thermal). Such thermal states have a characteristic function 

XT.(„,.)(0 = e-^'(^"°^'+^"'"') (2.53) 

for an average number of photons n and a squeezing parameter s (squeezing applied on x quadrature). Given 
similarly squeezed thermal states with the same average number of photons for both the external modes, the 
multiplying factor in eq. ( I2.52l i will be 

g - 2n^ e~^° sin(v3)^ Zg g - 2 sin{4>f w% (2.54) 

As explained before, the main difference with the examples of Gaussian mixtures in the previous sections 
is that this "smearing factor" will tend to become constant as the measurement gain goes to 1. 
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2.6 Teleportation fidelity with characteristic functions 

We have adopted the definition of fideUty proposed in ref . ||39l . In the characteristic function representation 
this is equal to eq. ( 11.71b . This is the scalar product, or the trace of the product of density operators, for the 
input and output states. Let us take, for simplicity of exposition, the output state we have derived in eq. ( 12.261 ) 
(for a maximum gain ~ gp = 1) and insert it in eq. (|1.71t ; 

^=^'' J xU^b) X^n{-^B) XAb{~Cb\ "^s) (2.55) 

A summary analysis might be made of this expression and of its integrand. We have inside the integrand 
that Xi„(^B) Xiri(— Cb) = |Xin(CB)P- TWs quantity is a real, even, nonnegative characteristic function itself, 
corresponding to the square of the density operator y5f„ (see eqs. (11.52b and ll.51b ). Its trace integral over is 
equal to the purity of the input state. In eq. ( 12.551 ) it is multiplied by the characteristic function of the resource 
state XAB- When the resource state approximates an EPR state, its characteristic function becomes nearly 
constant with a value of 1 (which is not square integrable). The result is ideal teleportation with maximal 
fidelity: 1 for pure input states and the value of the purity for mixed input states. 

Given that the fidelity is real, positive and bounded from above, we remark that xab(±^|j; ±^b) has to 
be bounded and must have an imaginary part that is odd or 0. Its real part must be even and positive (in the 
sense that an application is positive). 

The fidelity in eq. (12.55b is quadratic in the input state and linear in the resource state. The fidelities 
corresponding to the output states derived in sections l272l 12 . 3 1 and 12.51 share these two features; even though 
there are additional multiplying factors. 

One very simple conclusion to be drawn from the linearity of the teleportation output and of the fidelity 
with respect to the resource state is that we need only calculate the fidelity of teleportation for pure resource 
states, as the fidelity for a simple mixture of resource states comes about trivially. Given a mixture of pure 
resource states | ijji )ab, and given the characteristic function's linearity with respect to the density operator 
(see eq. (11.54b ). we have that 

PAB = '^Pi\lpt )aB ( \aB 

i 

Xab{Ca; Cb) = X! Xi>, ,ab{(,a; Cb) (2.56) 

i 

where Xi/>i , ab are the characteristic functions of the pure states | ipi )ab- 

Given the linearity of the fidelity (eq. (12.55b ) with respect to the above characteristic function; it follows 
that the fidelity of teleportation using this mixed resource state is simply the weighted average of the fidelities 
Ti resulting from the use of the pure resource states singly as teleportation resources; 

i 

Therefore, we need only concern ourselves with pure resource states, at least as regards discrete mixtures 
of resources. These discrete mixtures appear whenever one considers conditional preparation procedures for 
resources; such as photon subtraction and addition for two-mode squeezed Gaussian states, which involve 
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measurements. 

The fideUty of teleportation is quadratic with respect to the input state. For an input state that is a mixture 
of other (pure) input states with characteristic functions Xin,j{£.in), we have 



which is rather complicated. The integrals summed in eq. ( |2.59t are traces of the input pure state "i" with the 
output state resulting from the teleportation "f . The integrals are not (individually) bona fide teleportation 
fidelities; save for the case where i = j, but must sum to a teleportation fidelity that is bounded. The 
summation over the integrals is not equivalent to the average of fidelities in eq. ( 12.581 Even in the case of 
a near-ideal resource; eq. ( |2.59t includes terms that are the casual overlap of the components of the input 
mixture, and the (nearly) perfect transcriptions resulting from teleportation. This can be seen as a limitation 
in the concept of fidelity we have adopted. 

If we consider the mixture above described as the ensemble of likely (one qudit) inputs of a quantum 
teleportation channel for a CV 1681 or hybrid quantum computation scheme 1291 ; we can count on the 
likely inputs for the channel to be far away from each other as regards scalar product: nearly orthogonal 
or orthogonal, so as to be easily distinguishable by measurement apparatus even under noisy conditions. In 
that case, and for a good enough teleportation resource, we can consider eq. ( 12.591 ) to include mainly the 
contribution from the i — j traces (bona fide teleportation fidelities) that sum to a weighted average fidelity. 




(2.58) 



The fidelity of teleportation in eq. ( |2.55| for the above mixture is given by 




(2.59) 



Chapter 3 



Teleportation with degaussified, 
squeezed Fock and squeezed Bell-like 
resources 

In the previous chapter we introduced a formalism for the study of teleportation in CV based on the Wigner 
characteristic function representations of quantum states. The expression for the output state derived for the 
ideal CV protocol (section 12.21 ) is general for any combination of teleportation resource and input having 
Wigner and Wigner characteristic functions. 

The objective of the investigation of teleportation resources is to produce feasible teleportation schemes 
that improve the teleportation fidelity for the input states most likely to be used in CV quantum computa- 
tion 1681 120I or in hybrid schemes that use CV states for communication 1291 [3T1 . The original protocol |f33l| 
was developed for Gaussian two-mode squeezed states, and these will be used as the benchmark from which 
it is desirable to improve output state fidelity, for a fixed squeezing parameter r. The parameter r assumes 
the meaning of indicator of technological capability in the preparation of resources. 

An initial analysis of non-Gaussian resources can begin with the class of "degaussified" two-mode squeezed 
states that have been shown to produce an improvement in the teleportation fidelity with respect to Gaussian 
states of both a comparable co variance matrix and similar squeezing parameter r lISTl l22l l23l [3l . 

"Degaussification"Q is the production of a non-Gaussian state of the radiation field by conditional photon 
subtraction (addition) on an initially Gaussian state. This procedure is both performed and verified by 
appropriate single-photon measurements. The "simulation" of the conditional measurement and the use of 
the subsequent mixtures of degaussified and Gaussian resources has been done |§ in refs. l2ni23l [3l. Within 
the characteristic function representation, we have shown (see section l276l l that the use of any such mixture 
as a resource gives results for fidelity that can be trivially derived from the results for pure state resources. 

We have found useful for the search of improvements in teleportation fidelity to formulate classes of 
resource states that encompass the largest possible categories of Gaussian, degaussified or simply non- 
Gaussian states. These generalizing classes of resources are by definition non-Gaussian. They are necessarily 

' Not to be confused with the demagnetization procedure 

-In representations different from the characteristic functions'. 
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superpositions of the "special case" resources or involve additional unitary transformations on such states; 
moreover, they reduce to the special case states for particular choices of the superposition and unitary trans- 
formation parameters. These same parameters can be arbitrarily chosen, allowing us to sculpt teleportation 
resources; even to optimize the resource for maximal fidelity of teleportation of a given input state. A further 
and related advantage lies in the ability to compare optimal resources with special case resource states for 
which preparation strategies have been devised or experimentally tested; having been proposed in previous 
work for CV quantum information resources. 

In this chapter, we will first introduce the non-Gaussian resources and inputs used throughout the chapter, 
in the characteristic function representation. We will then study the fidelity in the teleportation of selected 
pure input states; using photon-subtracted and photon-added squeezed states, two-mode squeezed Gaussian 
states, two-mode squeezed Fock states; and as implied before, a class of states that includes the aforemen- 
tioned resource states as special cases. These we have named the squeezed Bell-like states in previous 
work 1241; general superpositions of the two-mode, Gaussian squeezed vacuum state and the two-mode 
squeezed Fock state of number 1, 1, which is non-Gaussian. For the squeezed Bell-like resource states, we 
have calculated an optimal fidelity (with respect to the superposition parameters) for each of the selected 
input states ll24l and compared this optimal fidelity with the (smaller) fidelities obtained for the resources 
mentioned above. 

The optimal squeezed Bell-like states and the "special case" resources are also studied and compared as 
to entanglement, non-Gaussian character and affinity with Gaussian states; with the purpose of establishing 
the characteristics of a non-Gaussian state that better correlate with teleportation fidelity. The parameters for 
comparison are the von Neumann entropy iTTSl ; a non-Gaussianity measure 11791 . and a two-mode squeezed 
vacuum "affinity" measure devised for non-Gaussian teleportation resources 1241 . 

3.1 Resources and inputs: state vectors and characteristic functions 

In this section, we will introduce the two-mode non-Gaussian states that are the object of our study as 
resources for CV teleportation in this chapter 

3.1.1 Squeezed Fock states 

The two-mode squeezed Fock state is prepared by the performance of a two-mode squeezing operation on a 
CV system that is in a (separable) two-mode Fock state. The state vector for this wholly non-Gaussian state 
is given by, 

IC ; TOA , TO_B ) = Sab{0 \ niA, mB )ab (3.1) 

where Sab (C) is the two-mode squeezing operator of eq. (11.411 1; with the separable state vector | , ) 12 
I ruA )a® I )b being the two-mode Fock state. 

The two-mode squeezed vacuum state Sab (C) I C ; 0; ) and the two-mode squeezed Fock state Sab (C) I C ; 1 
are special instances of this class of states. This last state will be referred to as two-mode squeezed Fock state 
throughout the chapter. 

The calculation of the characteristic function for the two-mode squeezed Fock state is straightforward. 
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given that it is pure state; hence it's characteristic function is given by eq. ( 11.55b . Thus 

X^SNiU;^B) = (1, 1| S\siO Da{U) Db{^b) Sab{0\ 1, 1) (3.2) 

Recall the Bogoliubov transformation effected by the unitary two-mode squeezing operator on the dis- 
placement operators (see eq. (|1.45t ) and the expression for the matrix element of the displacement operator 
in the Fock basis 11461 ; 

{m\m)\n) = (^y^'r^"e-il«l'4'"-")(|eP) (3.3) 

where il™""^ is the associated Laguerre polynomial lISOl . 

Given the above relations, the characteristic function for the state | C ; 1 , 1) is easily seen to be 

^^Ua, ) = ■/ ^ (I «^ 1^ + I 1^) (1 - I ^'^ n(l - I ^B ?) (3.4) 

where the variables and are related to the variables and via the Bogoliubov transformation 
described in eq. (11.461 1 for the displacement operator arguments. 

3.1.2 Degaussified resource states 

The degaussified resource states are generated by photon subtraction (or addition) to each of the modes A 
and i? of a two-mode squeezed Gaussian state. Although the photon subtraction (addition) procedure is 
conditional in most experimental setups, it will be assumed here that it has been successfully performed and 
verified, obtaining a pure resource state. Thus, the photon subtracted (added) states have the state vectors 

, m^-^ ■ C) -AAj-) a™- a^- Sab{0 |0, 0)^s (3.5) 
|m(+) , m^+^ ; C) =U^l^ a^^^ a^"- Sab{C) I , 0)^b (3.6) 

where are normalization constants specific to each state vector. Henceforth, and in keeping with the 
proposals for the experimental generation of these states ||7][8]|9][T0][lT], we will restrict our analysis to the 
case where niA = r?7,B = 1; to the single-photon subtracted (added) states. 

To easily calculate the normalization constants and the respective characteristic functions, and to put the 
degaussified states in a proper perspective with respect to the other resources studied in this chapter; we 
write the states described in eqs. ( I3.5l l and ( I3.6l l, taking into account the transformation that results from the 
application of the two-mode squeezing operator on the annihilation and creation operators (see eq. ( 11.44b ). 

1^"^ =AAe** 5^5(0 (-10,0)^5 + e'nanh(r)|l, (3.7) 
|l(+), C) -AAe-^-^^^sCC) (-tanh(r)lO, 0)ab + e'^ll, (3.8) 

where A/" = [ 1 + tanh^(r) is the normalization constant for both desgaussified states. Thus, degaussi- 
fied resources are shown to be superpositions of the two-mode squeezed vacuum and the two-mode squeezed 
Fock state. 

However, eqs. (13.7b and (13.8b substantially differ in the exchange of the hyperbolic coefficients for the 
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superposition, and in the character of the state for vanishing squeezing. Even if both states become separable 
for r = 0, the photon-added squeezed state reduces to the non-Gaussian two-mode Fock state; while the 
photon- subtracted squeezed state becomes the Gaussian the two-mode vacuum. 

The calculation of the characteristic functions for these degaussified resources makes use of the same 
relations (eqs. (11.45b . (11.461 1, (13.2b and (13.3b ) used in the derivation of the characteristic function of the two- 
mode squeezed Fock state (eq. (13.4b ). Recall that the degaussified resources are superpositions of two-mode 
squeezed Fock states (eqs. (13. 7b and (13. 8b ). The density matrices have the form 

p = iifonc,o,o)(c,o,oi + ixinc,i,i)(c,i,ii 

+ A'o*^i|C,l,l>(C,0,0| + K*i^o|C,0,0)(C,l,l| (3.9) 

with the factors A'o.i taking the appropriate values for photon-subtracted or added states. Terms in eq. ( 13.9b 
correspond to (Fock basis) matrix elements of the displacement operator (see eq. ( 13. 3b ) in the calculation of 
the characteristic function. The characteristic functions of the photon-subtracted and photon-added two-mode 
squeezed states are thus given by 

+ tanh2(r)(l-|e:in(l-KBp)} (3-10) 

xi'/iie^; ^b) =^2g-V2(l«:.l^ + iei,n {tanh2(r) - 2 tanh(r) Re[e- ^^^^i C^] 

+ (l-leAn(l-ieBl')} (3.11) 

where the relation of the variables and to the variables and is described by the Bogoliubov 
transformation of eq. ( 1 1 .46b for the displacement operator arguments. 

Comparing eq. (11. 58b with eqs. ( |3.4b , ( 13.10b , and ( 13.1 lb ; we see that the polynomial terms that define 
the non-Gaussian character of the state are always multiplied by a Gaussian factor equal to the two-mode 
squeezed state characteristic function of eq. ( 11.58b . 

Lastly; it is worth remarking that the (non-Gaussian) two-mode photon-subtracted squeezed state can be 
defined as the first-order truncation of the (Gaussian) two-mode squeezed state. Let us consider the two- 
mode squeezed vacuum given by | — 2 r ) = Sab{~ 2 r) | , 0)ab- RecalUng eq. ( 11.43b ; such a state can 
be written as 

\~2r)AB ^SAB{-r)SABi~r)\0,0)AB 

OO 

^SAB^r) cosh~i(r) ^ (tanh(-r))" \ n,n)AB (3.12) 

n=0 

Truncating the series of eq. (13.12b beyond n = 1 , we recover the photon-subtracted resource state 
(eq. (13.7b ), with (j) ^ n; that is r). Moreover, this state coincides with that of the photon- 

subtracted state introduced in ref. ||3l for the ideal case of a beam splitter with unity transmittance. 
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3.1.3 General non-Gaussian resources: Squeezed Bell-like states 

In order to unify the study of the properties, of the teleportation performance, and of the experimental methods 
for the generation of the above mentioned resource states, we take into consideration a class of states that have 
been named the squeezed Bell-like states in previous work [24 1. The state vector for the squeezed Bell-like 
state is given by 

I ■^)sBeU = Sab{0 ( cos{5) | , )ab + sin(5) 1 1 , 1 )ab) (3.13) 

which is a superposition of the two-mode squeezed vacuum and the two-mode squeezed Fock state. The su- 
perposition coefficients are parameterized as cos((5) and e' ^ sin((5), which is convenient as the free parameters 
i5 and 6 involved can be chosen arbitrarily; always having a normalized state vector 

The state in eq. ( 13.13b is two-mode squeezed. That it is Bell-like is apparent on inspection of eq. ( 11.641 ). 
For r ~ and 5 = 7r/4 the squeezed Bell-like state reduces to | , )ab + e'^ 1 1 , 1 )ab\ the Bell state for 
the two-dimensional Hilbert space spanned by Fock states | ) and | 1 ) ; and the ideal resource for discrete 
variables quantum teleportation lISTl of states belonging to that two-dimensional Hilbert space. 

The squeezed Bell-like state will have as special cases all of the two-mode states that we consider as 
teleportation resources in this chapter; for appropriate values of the parameters 5 and 6, and will interpolate 
between these resources. Taking ^ r, we obtain a Gaussian resource for i5 = 0, a squeezed Fock state for 
(5 = 7r/2, a photon-subtracted squeezed state for 5 = cos~^(A/'), 9 = tt and a photon-added squeezed state 
for S — cos^^(A/'tanh(r)), 6' = tt. The squeezed Bell-like state is inseparable for r = 0, except for the 
aforementioned trivial choices of 6 corresponding to squeezed vacuum and squeezed Fock state. 

The advantage of interpolating between known resource states allows for an unifying study of their char- 
acteristics and of their performance as entangled resources for quantum information in CV. This advantage 
is compounded by the possibility of choosing the parameters S and 9 in an arbitrary manner, thus sculpting 
the entangled resource. For instance, given an input state, and an analytical expression for teleportation 
fidelity; this fidelity can be optimized with respect to the superposition parameters, for a fixed squeezing r. 
Even if squeezed Bell-like states are considered nothing more than theoretical constructs for the study of the 
teleportation fidelity for a wide variety of resources; it will always be possible to pick a suitable "special case" 
to sculpt in an experimentally feasible manner so as to better approximate the optimal squeezed Bell-like state 
resource. 

The squeezed Bell-like state is, like the degaussified states, a superposition of two-mode squeezed states. 
The density matrix for the squeezed Bell-like state is of the general form described in eq. (13.9b . Calculating 
the characteristic function for this state involves repeating the procedure outlined above for the degaussified 
states; for a superposition with coefficients cos((5) and e*^ sin(5). The characteristic function of the squeezed 
Bell-like state reads 

XSBeiiiU; Cb) =e-i/^ (I + I') {cos2(5) + 2cos(<5) sin(<5) Re[e'''^:, 

+ sin2(5)(l - 1^:4^(1 - le^ni (3.14) 

where, as before; and are related to £^a and £,b by the Bogoliubov transformation of eq. ( 11.46b for the 
displacement operator arguments. 
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3.1.4 Selected input states 

The following single-mode input states have been considered for the study of teleportation with non-Gaussian 
resources, among a variety of Gaussian and non-Gaussian states that are likely inputs both in CV quantum 
computation and hybrid quantum computation: Coherent states; squeezed vacuum states; single-photon Fock 
states; squeezed single-photon Fock states and photon-added coherent states. 

The characteristic functions for most of these one-mode states are straightforward to derive, and can be 
found in the literature, for example in refs. 1471143)1 . 

For the coherent state | (3 ) the characteristic function reads 

Xco.(6„) =e-il«-l^ + 2»^-[f">^-] (3.15) 

The one-mode squeezed vacuum state | e ) = S{e) \ ) (where (e = e''^ s)) has a characteristic function 
given by 

X.gten) = e-^l«^"l' (3.16) 

where ^-^ = cosh(s) + f*„e"^ sinh(s), a Bogoliubov transformation of variables. 

The characteristic section for the squeezed Fock state S{£) \ 1 ) can be calculated easily given eq. ( I3.3l l 
and reads 

Xsgi=^ten)= e-^l«-"l' (1 - |e,'„n (3.17) 

with given by the above Bogoliubov transformation. 

The characteristic function for the Fock state 1 1 is given by the trivial limit s = of eq. (|3.17t . 
The (non-Gaussian) photon-added coherent state is prepared by the application of the creation operator 
to a coherent state; 

(1 + |/3|2)-i/2 5t|^) (3.18) 

where, for /? = eq. ( 13.181 ) reduces to the single-photon Fock state | 1 ) . The derivation of the characteristic 
function of the photon-added coherent state is; however, straightforward 

Xpac{^^n) = (1 + | | ^ ) - 1 1 1 + ^ M?-. ^ 1 (1 + „ |^^j2 ^ 2jIm[C„/3*]) (3.19) 

3.2 Teleportation with degaussified and squeezed Fock resources 

In this section, we compare the behavior of the teleportation fidelity for the input states described above; 
using as resources the two-mode squeezed Fock state (eq. (I3.1I) ). and both the degaussified states (eqs. (I3.7I) 
and dlSTl). 

We will assume ideal CV teleportation. Therefore, the output state characteristic function is given by 
eq. (12.261 1 (for = Sp = 1), and the teleportation fidelity is given by eq. ( 12.55b . The fidelity is analytically 
computable for the the input states and entangled resources considered here; the integral in eq. ( 12.551 1 can be 
exactly calculated in terms of finite sums of averages over Gaussian distributions. 

The fidelities of teleportation will be analyzed and compared for fixed and equal squeezing parameter 
r for all resources. The performance of the Gaussian two-mode squeezed state will also be displayed for 
reference purposes, as it is generally considered a benchmark 1391 l40l . This choice is made to compare the 
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performance of Gaussian and chosen non-Gaussian resources, given the possession of the same technological 
means for their generation. 

3.2.1 Teleportation of Gaussian input states 

We begin our analysis with the teleportation of Gaussian input states. In fig. (13.2.1b we plot the fidelity 
of teleportation (see eq. dl.Tll i) having as inputs the coherent state | /3 )m (Panel I), and the input squeezed 
vacuum state | jS ) ;„ (Panel II), for the resource states introduced above, with the exception of the squeezed 
Bell-like state. The phase of squeezing will be fixed as (/) = tt from now on; in keeping with the conventions 
established in the CV protocol outlined in chapter|2] 

For both inputs in fig. (13.2.1b ; the choice of the photon-subtracted squeezed state (eq. (13.10b ) as entangled 
resource is the most convenient one. It corresponds to the highest value of the fidelity T for any fixed value of 
the squeezing r in the realistic range [0, 1]. On the contrary, the choice of the squeezed Fock state (eq. ( 13.4b ) 
as entangled resource is the least convenient, yielding the poorest performance; even when compared to the 
Gaussian squeezed resource. Finally, regarding the use of the photon-added squeezed state (eq. ( |3.11b ) as 
entangled resource, it allows for a very modest improvement in the fidelity when compared to the Gaussian 
resource; and this only for a small interval of values around r = 1. For smaller squeezing, its performance is 
poorer than that of the Gaussian resource. 

The poor performance of the photon-added and two-mode squeezed Fock resources is not surprising for 
small squeezing parameters. The input in both cases is a Gaussian state; and both these resource states reduce 
to a two-mode Fock state for r = 0. In Panel I of fig. (13.2.1b we see that the overlap of a coherent input 
state with a mixture of displaced^ coherent states for classical teleportation is 0.5, as expected from previous 
work on "classical" teleportation ||39ll40l . The overlap of the coherent input state with a mixture of randomly 
displaced one-photon Fock states is 0.25. The probability for a given displaced element of the mixture is given 
by the probability of a result equal to the displacement; when the joint measurement of non-local quadratures 
of the A and in modes is performed (see subsection ll.4.2b . 

We must, finally, remark that the fideUty of teleportation is significantly higher for the coherent state input 
(Panel I,fig. ( 13.2.1b ) than for the squeezed vacuum input (Panel II,fig. (13.2.1b ). for all resources, for small r. 
This is due to the fact that the one-mode squeezing operation has an asymptotic limit s ^ oo in which 
all squeezed states will approximate a quadrature eigenstate. In the case of quasi-classical teleportation, for 
small r, the randomly displaced, nearly uniform mixture in mode B has little overlap with a squeezed input 
state; even for the moderate squeezing (s ~ 0.8) chosen. 

3.2.2 Teleportation of non-Gaussian input states 

Let us now consider the teleportation performance for non-Gaussian input states; the single-photon Fock 
state I 1 )in, the squeezed single-photon Fock state Sin{s) \ 1 and the photon-added coherent state (1 + 
I p'j-i/a I ^ ^ In fig (13. 2. 2b . we plot the fidelity of teleportation for two non-Gaussian input states: the 
single-photon Fock state (Panel I), and the photon-added coherent state (Panel II). 

In Panel I, we observe that both the degaussified resources lead to an improvement of the fidelity with 
respect to the Gaussian resource. The photon-subtracted squeezed resource again performs better than the 
photon-added resource; while the squeezed Fock state yields the poorest performance when compared to 

^Randomly displaced, given the separability of the "resource" 
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Figure 3.1: Fidelity of teleportation JF, as a function of the squeezing parameter r, with (/> = tt, for 
input coherent state |/3) (Panel I) and input squeezed vacuum state | e ) (Panel II). Comparison is given for 
different two-mode entangled resources; {a) squeezed state (full line); {h) squeezed Fock state (dashed line); 
(c) photon-subtracted squeezed state (double-dotted, dashed line); (d) photon-added squeezed state (dot- 
dashed line). In plot I the value of /3 is arbitrary. In plot II the squeezing parameter of the input state is fixed 
at s = 0.8. 
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Figure 3.2: Behavior of the fidelity of teleportation as a function of the squeezing parameter r, with 
(j) ~ TT, for two different non-Gaussian input states: The Fock state |1) (Panel I), and the photon-added 
coherent state (1 + |/3|2)-i/2at|/3) (Panel II). We compare the performances obtained by using different two- 
mode entangled Gaussian and non-Gaussian resources: (a) squeezed state (full line); (b) squeezed Fock state 
(dashed line); (c) photon-subtracted squeezed state (double-dotted, dashed line); (d) photon-added squeezed 
state (dot-dashed line). In Panel II the value of the coherent amplitude of the input photon-added coherent 
state is fixed at /3 = 0.3. 
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all the resources. From Panel II we see that once more the photon-subtracted resource yields the best 
performance for any fixed squeezing parameter; that the photon-added squeezed resource allows for a very 
modest improvement in the fidelity with respect to the squeezed Gaussian reference for small squeezing. 

Differently from the teleportation of Gaussian inputs (see fig. ( 13.2. It ), the photon-added resource is no 
worse than the Gaussian state for small squeezing r. That the two-mode squeezed Fock state is a poor 
resource even at a small squeezing seems unusual. Let us recall that the overlap in eq. ( 12.551 1 for classical or 
nearly classical teleportation (r — > 0) in this case is that of a non-Gaussian, pure input state with a mixture 
of randomly displaced Fock states; not with a single-photon Fock state. 

It is interesting to note that teleportation fidelities (for the same squeezing value r) are higher for all the 
resource states when the input is a photon-added coherent input (compare Panels I and II of fig. (13.2.2b ); 
higher than when the input is the Fock state that is the zero amplitude limit of the photon-added coherent 
state. The improvement is small, however, for the input and resource states we have considered. This is not 
unsurprising, given that the photon-added coherent state, at amplitude |/3p 3> 1 approximates a coherent 
state. The protocol for CV teleportation seems more suitable for the teleportation of coherent states than for 
the teleportation of Fock states, at least as far as our chosen definition of fidelity is concerned. 

In fig. ( 13.2.2b we compare the fidelity of teleportation T for the case of a squeezed Fock input state and 
different Gaussian and non-Gaussian resources. Comparing with Panels I and II of fig. (13.2.2b . we see that 
the qualitative behavior of teleportation is different from both previous examples, and more reminiscent of 
that for the Gaussian inputs (see fig. ( 13.2.1b ). This we explain by noting that the one-mode squeezed Fock 
state is (at s = 0.8) of a more Gaussian character than the previous examples of Fock state and photon-added 
coherent state at small amplitudes. 

li I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' 
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Figure 3.3: Behavior of the fidelity of teleportation as a function of the squeezing parameter r, with 
= TT, for the squeezed Fock input state S{s) \ 1 ), using different two-mode Gaussian and non-Gaussian 
entangled resources: (a) squeezed state (full line); (b) squeezed Fock state (dashed line); (c) photon- 
subtracted squeezed state (double-dotted, dashed line); (d) photon-added squeezed state (dot-dashed line). 
The value of squeezing is for the input state is fixed at s = 0.8. 

From all of the above investigations, we conclude that the photon-subtracted squeezed state (eq. ( 13.7b ) 
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is always to be preferred as entangled resource compared either to the Gaussian ones or to non-Gaussian 
states that are obtained by combining squeezing and photon pumping. The reason explaining this result will 
become clear in the next sections when we will discuss the general class of states that include as particular 
cases all the resources introduced so far; and attempt to single out some properties of resource states that are 
related to improved performance in quantum teleportation. 

3.3 Squeezed Bell-like resources for teleportation 

In this section we study the fidelity of CV teleportation of the squeezed Bell-like state of eq. ( |3.13l l. The 
class of squeezed Bell-like states includes as special cases all the resources studied so far, as well as the Bell 
states of the two-dimensional Hilbert space spanned by the vacuum and single-photon Fock states. Part of 
the analysis will consist in the optimization of the teleportation fidelity with respect to the free superposition 
parameters of squeezed Bell-like states. The optimal squeezed Bell-like state (for a given input state) can 
always point the way to an adequate (if not optimal) choice of technically feasible resources for teleportation; 
and to the further sculpturing of such resources. Even if some states in the class of Bell-like states remain out 
of the reach of experimental realization. Moreover, the study of the properties of a wide class of resources can 
be unified and interpolated by the use of squeezed Bell-like states, not only regarding teleportation fidelity, 
but other uses as well. 

The teleportation fidelity (eq. (12.55b ) is computable in an analytical manner for all the input states 
considered, using the squeezed Bell-like resource. It is an expUcit function T{r, (j>, S, 9) of the independent 
parameters r, cf), S and 9 that describe the squeezed Bell-like resource. 

We do not report here the explicit analytic expressions of the fidelities associated to the squeezed Bell-like 
resource and to each input state, because they are rather long and cumbersome. Some of them are reported in 
Appendix lAl 

An explicit analysis has estabUshed that all the fidelities are monotonically increasing functions of the 
squeezing parameter r at maximally fixed phase (j> = tt. In the following we assume that (p = tt and 
6* = 0. It can be checked that non vanishing values of 9 do not lead to an improvement of the fidelity in CV 
teleportation. 

3.3.1 Teleportation with a squeezed Bell state 

At finite squeezing r and for S = j and 9 — 0, state (eq. (13.131 )) reduces to a squeezed Bell state. This 
resource is, as we have remarked in subsection 13 . 1 .31 intrinsicallv nonclassical; as well as having as a r = 
limit the Bell state for the two-dimensional Hilbert space spanned by the vacuum and the single-photon Fock 
state. 

We may assess analytically the performance of such an entangled resource as far as teleportation is 
concerned. In fig. ( 13.3. Il l we show the behavior of the fidelity as a function of the squeezing parameter 
r, with (j> = TT, S = J and 9 = for the five different input states considered in the previous Section. It 
is straightforward to observe that the squeezed Bell state; used as an entangled resource, leads to a relevant 
improvement of the performance, when compared to all the other Gaussian and non-Gaussian resources that 
we have investigated in the previous section. 

The squeezed Bell state will always have a better performance for Gaussian inputs than non-Gaussian 
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inputs. The "best" non-Gaussian input, the small amplitude = 0.09) photon-added coherent state, 

shows a noticeable improvement in fidelity with respect to both Fock and squeezed Fock states. This suggests 
that the CV teleportation protocol will always have a better performance for Gaussian (or approximatively 
Gaussian) inputs. 
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Figure 3.4: Behavior of the fidelity of teleportation JF(r, S, 9) associated to the squeezed Bell resource 
with (j) ~ -K, 5 = 9 = Q, plotted as a function of the squeezing parameter r for the following input 
states: (a) coherent state (full line); (b) squeezed state \s) = S{s)\Q), with s ~ 0.8 (dotted line); (c) Fock 
state |1) (dashed line); (d) photon-added coherent state (1 + \l3\'^)~'^f'^a^\P), with /3 = 0.3 (dot-dashed fine); 
(e) squeezed Fock state \s) — S{s)\l), with s = 0.8 (double-dotted, dashed fine). 

3.3.2 Optimized squeezed Bell-like resources for teleportation 

We proceed to maximize, for every input state, the fidelity T{r, tt, S, 0) over the Bell-superposition angle S. 
At a fixed squeezing r = f, we define the optimized fidelity as 

^opt{f) = max J-{f, TT, (5, 0) . (3.20) 

For a coherent state input the maximization of JF( r, tt, (5, ), at fixed r, leads to the following determina- 
tion for the optimal Bell-superposition angle Smax', 

•^IJ^L = ^ arctan(l + (3.21) 
For a single-photon Fock state input, the optimal angle is given by 

= 2 ^'^ct^^ \ 3(g2r_i)2 j (3.22) 

We report, in fig. (13.3.2b the behavior of the optimized fidehties Toptir) as functions of r for all the 
input states considered in this chapter. In all: coherent state, squeezed vacuum, single-photon Fock state. 
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Figure 3.5; Plot of the fidelity Toptir), using optimal squeezed Bell-like resources, as a function of r for the 
following input states: (a) coherent state (full line); (b) squeezed vacuum | s), with s = 0.8 (dotted line); 
(c) single-photon Fock state 1 1 ) (dashed line); (d) photon-added coherent state (1 + |/?p)^^/'^d^|/3), with 
/? = 0.3 (dot-dashed line); (e) squeezed Fock state |1, s ), with s = 0.8 (double-dotted, dashed line). 

A relevant improvement of the fidelity is observed in all cases, even at vanishing squeezing, due to the 
persistent nonclassical character of the optimal, squeezed Bell-like entangled resource; even in the limit 

The Gaussian states, and the photon-added coherent state show higher teleportation fidelities than the 
wholly non-Gaussian states for all r values. This result; together with similar results in subsections 13.2.21 
and l3.3.ll stronglv suggest that some local operations; a suitable displacement, for instance on a non-Gaussian 
input state previous to CV teleportation might improve the fidelity of teleportation, for any resource used. 

To quantify the increase in the probability of success for teleportation when we use the optimal squeezed 
Bell-like state for a resource; we define a relative fidelity coefficient AjF(r); relative with respect to the 
fidelity Tref {r, tt) achieved using with one of the "special case" resources as a benchmark . A natural 
definition of such a relative quantity is 

Tref [r, tt) 

In fig. ( I3.3.2I I, we plot the relative fidelity A J-'{r) as a function of the squeezing parameter r for two 
different choices of benchmark resources that we have deemed the most interesting, given the high values 
of the fidelities obtained when using them. In Panel I, the benchmark resource is the Gaussian two-mode 
squeezed vacuum; in Panel II the benchmark resource is the non-Gaussian two-mode photon-subtracted 
squeezed state. 

In Panel I of fig. ( I3.3.2l i we see that at fixed squeezing, the use of optimized squeezed Bell-like resources 
leads to a strong percent enhancement of the teleportation fidelity (up to more than 50%) with respect to that 



50 



Chapter 3 Teleportation with degaussified, squeezed Fock and squeezed Bell-like resources 




Figure 3.6; Behavior of the relative fidelity AjF(r) using an optimal squeezed Bell-like resource, as a function 
of r for the following input states: (a) coherent state (full line); (b) squeezed state \s), with s = 0.8 (dotted 
line); (c) single-photon Fock state |1) (dashed line); (d) photon-added coherent state (1 + \j3\'^)^^^'^a)\j3), 
with (3 = 0.3 (dot-dashed line); (e) squeezed Fock state |1, s), with s = 0.8 (double-dotted, dashed line). In 
Panel I the benchmark resource is the two-mode squeezed vacuum; in Panel II the benchmark resource is the 
two-mode photon-subtracted squeezed state. 
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attainable exploiting the standard two-mode squeezed vacuum. The greatest relative improvements, and the 
maxima in relative fidelity, occur in the teleportation of the non-Gaussian inputs, particularly with the Fock 
state and the photon-added coherent state (for which the relative improvement is not as great). Obviously, in 
the asymptotic limit of very large squeezing, the two resources converge to unity teleportation efficiency, as 
the resources themselves converge to EPR states. 

Panel II shows that the use of the optimized squeezed Bell-like resource leads to a significant advantage 
with respect to the use of the photon-subtracted squeezed state resource for low values (up to r ~ 0.5) of 
the squeezing. The different curves corresponding to the different input states, exhibit the same qualitative 
behavior; however, the greater relative improvements are those for the fidelity of teleportation of non- 
Gaussian inputs, particularly the Fock state and the photon-added coherent state (again, with a lesser relative 
improvement than the Fock state). Starting from large values, AT{r) decreases monotonically and vanishes 
at different points in the interval [0.5 < r < 0.9]. After vanishing, the relative fidelity exhibits peaks 
at different intermediate values of the squeezing, before vanishing asymptotically for large values of r. 
Unsurprisingly; for these minima values of r = f such that AjF(f) = 0, the optimized Bell-like state and the 
photon-subtracted squeezed state are identical. 

3.4 A comparison of entanglement, non-Gaussian character and Gaus- 
sian affinity of resource states 

In this section, we analyze the characteristics of the resource states studied in this chapter, to determine which 
properties of the resources have an influence on teleportation fidelity. Particularly, we study the entanglement 
and non-Gaussian character of optimized Bell-like states and compare them with those of photon- subtracted 
and photon-added two mode squeezed states. This comparison is to be particularly illuminating, given the 
difference in teleportation fidelity obtained using these degaussified states as resources; and given the very 
different characteristics (obvious for r — ^ 0) of these two types of resources that seem otherwise similar 

First, we study entanglement, as given for pure states by the von Neumann Entropy. Then, we study non- 
Gaussian character as given by a measure of the distance between the resource under study and a reference 
Gaussian state having the same covariance matrix and first order averages as the resource. Thirdly, we 
define and study Gaussian affinity as given by the overlap of the resource and a two-mode squeezed vacuum, 
maximized over the squeezing parameter of the latter. 

We analyze these measures in conjunction for the resources above mentioned, paying special attention to 
the behavior of these measures for optimized Bell-like states. 

3.4.1 The von Neumann entropy of non-Gaussian resources 

The bipartite entanglement; the Schmidt rank of the pure non-Gaussian states chosen for study in this chapter 
can be quantified in an unique manner using the partial Von Neumann entropy (entropy of entanglement) 
EyN. Given the density operator for the two-mode state p, the von Neumann entropy is equal to the 
Shannon ifTSl entropy of the partial traces: 

E.N = -Tr (p^^(^) log,(p^^(^))) = -Tr log,(p^^(^))) (3.24) 
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where the partial traces are indicated by /5 = Tr^(/5), and by p-^'"^^' = Ti-b{p) 

To simplify calculations involving powers of the density operator of the aforementioned resource states 
(see section 13.11 ), the series expansion for the logarithmic function in eq. ( |3.24| i has been truncated to the 
first-order 

For the two-mode squeezed Fock and the degaussified states, the von Neumann entropy depends only on 
the modulus r of the squeezing parameter It is plotted in fig. (13.4.11 ) and compared to that of the Gaussian 
two-mode squeezed state. 
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Figure 3.7: Behavior of the von Neumann entropy E^j^ for two-mode squeezed Fock, two-mode photon- 
added and two-mode photon subtracted squeezed states, as a function of the modulus squeezed parameter 
r. The upper curve (dot-dashed line) corresponds to the squeezed Fock state |r ; 1 , 1 ); the intermediate 
curve (dashed line) corresponds equivalently to the photon-subtracted squeezed state r) and to 

the photon-added , 1'^+' ; r) squeezed state . The lower curve corresponds to the Gaussian two-mode 
squeezed state | r ) . 

At a given squeezing, all the non-Gaussian states show an entanglement larger than that of the Gaussian 
two-mode squeezed vacuum. In the range of experimentally realistic values < r < 1 of the squeezing, 
the two-mode squeezed Fock state is the most entangled state. Moreover, the photon-added and the photon- 
subtracted two-mode squeezed states exhibit exactly the same amount of entanglement at any r. 

However, the results in subsections 13.2. II and 13.2.21 indicate that, for a given squeezing r, a higher von 
Neumann Entropy (and greater entanglement) does not correspond to a higher teleportation fidelity. The 
two-mode squeezed Fock state obtains the smallest teleportation fidelities of any resource while having the 
greatest von Neumann entropies. And the photon-added two-mode squeezed state obtains lower teleportation 
fidelities than the photon-subtracted state, even if they are indistinguishable by their von Neumann entropy. 

The von Neumann entropy of the general squeezed Bell-Uke states is plotted in fig. ( 13.4. II ). In panel I we 
plot i?t,7v as a function of r and 6. In panel II, we can observe how the regular, oscillating behavior of the 
entropy for the Bell-like state [r = 0) becomes gradually deformed by the optical pumping (r 7^ 0), leading 
to a peculiar pattern of correlation properties for the squeezed Bell-like state. 

Note that for S ^ Q and S = 7r/2, we have the von Neumann entropy of the Gaussian state and the 
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squeezed Fock state, respectively in Panel II of fig. ( 13.4.1b . 

The squeezed Bell-like states have maxima of von Neumann Entropy for r ^ near 7r/4 and 3tt/A. 
However, it suffices to evaluate eqs. (|3.21t and ( 13.221 ) for r = (recall that the Bell-like state is entangled 
even if it is not squeezed) to realize that, while the optimal fidelity for the teleportation of Fock state inputs 
is obtained for Smax{r = 0) = 7r/4 (the Bell state of discrete variables settings), this is not the case for 
the coherent state inputs, with 6rnax{r = 0) = 0.554. Even in the case of Bell state resources for CV 
teleportation, the greatest von Neumann entropy does not maximum determine fideUties for certain classes of 
states. 

3.4.2 Optimized squeezed Bell-like resources: A comparison of von Neumann en- 
tropies 

In fig. (|3.4.2| i. we show the behavior of the von Neumann entropy E'^jv for two different squeezed Bell- 
like resources; the first one optimized for the teleportation of an input coherent state, with Smlixir) (see 
eq. ( I3.21l i): the second one optimized for the teleportation of a single-photon Fock state, with Smdxir) (see 
eq. (|3.22t ). This behavior is compared with that of the von Neumann entropy of both the degaussified states. 
The intersections between the curves in the figure correspond to the values os squeezing f for which the 
squeezed Bell-like state reduces to a photon-subtracted or to a photon-added squeezed state. 

It can be noticed in the range < r < f; in which the fidelity of teleportation using optimized Bell-like 
resources is substantially higher than that for the photon-subtracted squeezed resource (see fig. (I3.3.2l i. Panel 
II), that the entanglement (quantified by the von Neumann Entropy) of the squeezed Bell-like state is always 
larger than that of the photon-subtracted (and photon-added) squeezed states. Therefore, a partial explana- 
tion of the better performance of squeezed Bell-like resources lies in their higher degree of entanglement 
compared to other non-Gaussian resources at lower levels of squeezing. However, from the graphs it can be 
seen that for r > f the entanglement of photon-subtracted and/or added resources is larger; nevertheless, the 
fidelity of teleportation is below that associated to optimized squeezed Bell-like resources. Much lower, in 
the case of the photon-added resource, we remark. 

Therefore, entanglement is not the only parameter necessary to determine the teleportation performance 
of different non-Gaussian resources. The example chosen here validates this conclusion specially well; 
nonetheless because the squeezed Bell-like state is a generalization of degaussified resources, that reduces to 
a degaussified resource for certain values of the superposition angle S. 

3.4.3 The non-Gaussianity: A character measure 

We have concluded that entanglement alone does not suffice to determine the teleportation fidelity for a non- 
Gaussian resource. We have seen that the two-mode squeezed Fock state, the wholly non-Gaussian resource 
having the highest entanglement of the resources studied in this work, performs poorly as a teleportation 
resource. We have seen that the degaussified states, resources of a different character; but equal entanglement, 
perform the teleportation of the same input with very different fidelities. We have seen that the optimal choices 
of squeezed Bell-like resources for teleportation (at a fixed squeezing r) do not always correspond to maximal 
entanglement. We have seen notable differences in the teleportation fidelities of Gaussian and non-Gaussian 
inputs, using the same resource. 
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Figure 3.8: The von Neumann entropy E^n for the squeezed Bell-like state as a function of its defining 
parameters r and S. Panel I displays the three-dimensional plot of E^jy. Panel II displays two-dimensional 
projections at fixed squeezing strength r. Curves from bottom to top correspond to the different sections of 
EyN as functions of (5 for r = 0, 0.2, 0.4, 0.6, 0.8, 1. 
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Figure 3.9: The von Neumann entropy E.^n for the optimized squeezed Bell-Hke state, as a function of 

r. Dashed Une; Optimized resource for a coherent state input, 5 ~ slnax'-, Long dashed line: Optimized 

( p) 

resource for a Fock state input, S = 5„iax- The von Neumann entropy of the degaussified states is reported 
for comparison purposes (dot-dashed line). 



It is natural, as we study non-Gaussian resources and observe improvements in teleportation fidelity over 
Gaussian resources, to look for a quantification of the non-Gaussian character of the different resources we 
have considered, and in general of the optimized squeezed Bell-like resources in order to compare their 
performance in teleportation with respect to both this quantity and the entanglement in conjunction. The task 
is to define a reasonable measure of non-Gaussianity that is endowed with some nontrivial operative meaning. 

Recently, inspired by work |&| on the extremal nature of Gaussian states at fixed covariance matrix; a 
measure of non-Gaussianity has been introduced in terms of the Hilbert-Schmidt distance between a given 
non-Gaussian state and a reference Gaussian state with the same covariance matrix 1791 . Given a generic 
state with density operator p; its non-Gaussian character can be quantified through the distance dnc between 
p and the reference Gaussian state defined according to the following relation: 

_ Tt{{p- PG?) _ Tr(/3^) + Tr{pl) - 2T,{ppg) 

2Tr(/52) ~ 2Tr(p2) ' ^^'^^^ 

where the Gaussian state pc is completely determined by the same covariance matrix and the same first order 
average values of the quadrature operators associated to state p. 

Using this definition, in fig. ( 13.4. 3l l we report the behavior of the non-Gaussianity dnG for the squeezed 
Bell-like state. 

The quantity dno depends only on the superposition angle 5 of the squeezed Bell-like state (see Panel I), 
as the non-Gaussianity of the state cannot change under squeezing operations; the squeezing induced on the 
state measured would translate into squeezing of the reference Gaussian state. For 5 in the interval [0, tt], dnc 
attains its maximum at 5 = ^: At that point, the squeezed Bell-like state reduces to a two-mode squeezed 
Fock state. It is expected for a (two-mode squeezed) Fock state to be more non-Gaussian than a (two-mode 
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Figure 3.10: Non-Gaussianity measure dnc for the squeezed Bell-like state. Panel I shows dnc for the 
squeezed Bell-like state as a function of 8, for arbitrary r. Panel II shows dno for the squeezed Bell-like 



(dashed 



state as a function of r, for (5 fixed at the optimized values for coherent state inputs i5 
line); and b = 6max for Fock state inputs (long dashed line). For comparison, the values of dnC for the 
photon-subtracted resource (dot-dashed Une) and photon-added resource (double dotted, dashed line) are 
also reported. 
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squeezed) superposition of the vacuum and of a Fock state. 

In Panel II, we plot the behavior of dnc for the squeezed Bell-like resources optimized for the teleporta- 
tion of a coherent state input and a single-photon Fock state input, respectively, with S = Smax and 5 = Smlx- 
For comparison, we plot as well the non-Gaussianity dnc for the photon-added and the photon-subtracted 
squeezed states. The intersection points occur once again at the points f where the optimized squeezed Bell- 
like states reduce to the photon-subtracted squeezed states. For r in the range [0, f], the optimized squeezed 
Bell-like resources are not only more entangled; they are more non-Gaussian than the photon-subtracted 
squeezed States. We note that for lim^-foo fJrna^: = ^vair^oo Smax ~ 1. Thus, for very large squeezing the 
two optimized squeezed Bell-like resources tend to the state S'y!iB(—r) {cos^|0, Q) ab + sin ^ | 1 , 1)ab}, 
exhibiting a dominating Gaussian component. On the other hand, for large r, the squeezed photon-added and 
photon-subtracted squeezed states asymptotically tend to a squeezed Bell state (corresponding to Smax = f ), 
which has balanced Gaussian and non-Gaussian contributions. While for smaller values of r, they have very 
different values of non-Gaussianity; which is not surprising if we recall the enormous difference in character 
of these states for r — > 0. 

It is also remarkable that, for r in the range [0, f], that the optimal squeezed Bell-like resources for 
teleportation of the non-Gaussian Fock state {5max) have a much higher non-Gaussianity than those optimal 
squeezed Bell-like resources for teleportation of the Gaussian coherent state (Smax)- This might be under- 
standable for completely "classical" teleportation with an squeezed resource of the type studied so far, where 
fidehty is just the overlap of the input state in with a mixture of states constituted of different instances of 
random displacements of mode B. But this is not the case here; the squeezed Bell-like resources are always 
entangled even for r = 0, with the exception of some trivial choices for S. 



3.4.4 The Gaussian resource affinity 

We have developed a non-Gaussianity measure for entangled resources that is in effect a distance between the 
resource state we are studying and a Gaussian state of the same covariance matrix. But this Gaussian reference 
is different for every resource studied, and the interest of such a measure lies in making comparisons of the 
non-Gaussianity for different non-Gaussian resources. Moreover, we do not know that the Gaussian state 
constructed for every particular measure is a resource for teleportation, for it is not described as such. 

Observing that the squeezed Bell-like states and the photon-added and photon-subtracted squeezed states 
are all obtained through a degaussification protocol from a pure squeezed state, one could modify the defini- 
tion of eq. ( 13.251 1 by taking the two-mode squeezed vacuum | C' )ab {C = ) as the universal reference 
Gaussian state rhoa- This choice makes it possible to compare our non-Gaussian resources to an unique 
reference that is also a well studied teleportation resource. This will be particularly important if we are to 
analyze measures of non-Gaussianity, teleportation fidelities and von Neumann entropies together 

Adopting this modified definition, and observing that the non-Gaussian states to be compared and the 
reference Gaussian state are all pure, eq. ( 13.251 ) reduces to dnG = minr'^^'jl — Tr[ppG]], where the 
minimization is constrained to run over the squeezing parameters Q' of the reference twin-beam. However, 
it turns out that this modified definition provides results and information qualitatively analogous to those 
obtained by applying the original definition. 

There is still one, and apparently opposite property that plays a crucial role in the sculpturing of an 
optimized non-Gaussian entangled resource. From figs. ( 13.4.21 ) and ( 13.4.3b we see that at sufficiently large 
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squeezing the photon-added and photon-subtracted squeezed resources have entanglement comparable to 
that of the optimized squeezed Bell-like states and, moreover, possess stronger non-Gaussianity. Yet they are 
not able to perform better than the optimized Bell-like resources. This fact can be understood as follows, 
leading to a definition of squeezed vacuum affinity: It is well known that the Gaussian two-mode squeezed 
state in the limit of infinite squeezing realizes the EPR state which is the equivalent of the Bell (eq. (11.64b ) 
state for two-dimensional systems. These two ideal resources, respectively in the CV and discrete variables 
case, allow perfect quantum teleportation with maximal fidelity. Therefore, we argue that, even when 
exhibiting enhanced properties of non-Gaussianity and entanglement, any efficient resource for CV quantum 
information tasks should enjoy a further property, i.e. to resemble the form of a two-mode squeezed vacuum 
as much as is possible, in the large r limit. The squeezed vacuum affinity can be quantified by the following 
maximized overlap: 

g = ma.x\AB{-s\iljres{r))AB\'^ (3.26) 

s 

where | — s)ab is a two-mode squeezed vacuum with real squeezing parameter — s, and |VVes(7'))AB 
is any entangled two-mode resource that depends on the squeezing parameter r. This definition applies 
straightforwardly to the photon-added and photon-subtracted squeezed resources, and as well to the squeezed 
Bell-like resources whenever they are optimized with respect to a particular input state. The maximization 
over s is imposed in order to determine, at fixed r, the two-mode squeezed vacuum that is most affine to the 
non-Gaussian resource being considered. 

In fig. ( 13.4.41 ) we study the behavior of the overlap Q as a function of the squeezing r for different non- 
Gaussian entangled resources. 
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Figure 3.11: Maximized overlap G between a two-mode squeezed vacuum and different non-Gaussian 
entangled resources as a function of r: squeezed Bell-like state with delta fixed at the optimized value 
S = Smax for coherent state inputs (dashed line); the same state with delta fixed at the optimized value 
S ~ Smax for Fock state inputs (long dashed line). For comparison purposes the maximized overlap with 
the photon-subtracted squeezed state (dot-dashed line); photon-added squeezed state (double-dotted, dashed 
line); and the single-photon squeezed Fock state (dotted line) are displayed as well. 
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Comparing fig. ( 13.4.31 ) and fig. ( 13.4.4b we can see that the behaviors of the non-Gaussianity and squeezed 
vacuum affinity seem to be complementary. Take the curve for any one state and observe growth and 
asymptotic behaviors. For example the photon added resource: non-Gaussianity approaches its asymptotic 
value from above; it is obviously a convex function; two-mode vacuum affinity approaches its asymptotic 
value from below, and it is a concave function. 

Observe (in fig. (I3.4.4l i) that the photon-subtracted and photon-added resources that have equal entangle- 
ment for every r are perfectly distinguishable by squeezed vacuum affinity. 

In fig. (|3.4.4t . the greatest affinity Q is always achieved at large values of the squeezing parameter by the 
optimized squeezed Bell-like resources, while the lowest, constant affinity is always exhibited by the squeezed 
Fock states. Observe additionally that the Squeezed Bell-like resource optimized for coherent states (Sm\x) 
has a higher affinity than the resource optimized for Fock states {{Smax)), specially for the smaller values of 



In conclusion, optimized squeezed Bell-like resources are the ones that in all squeezing regimes are closest 
to the simultaneous maximization of entanglement, non-Gaussianity, and affinity to the two-mode squeezed 
vacuum. The optimized interplay of these three properties explains the ability of squeezed Bell-like states 
to yield better performances, when used as resources for CV quantum teleportation, in comparison both to 
Gaussian resources at finite squeezing and to the standard degaussified resources such as the photon-added 
and the photon-subtracted squeezed states. 
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While two-mode Gaussian squeezed states are currently produced in the laboratory, the experimental gener- 
ation of non-Gaussian states in quantum optics is still a complex task, as it requires the availability of large 
nonlinearities and/or the arrangement of proper apparatus for conditional measurements. Nevertheless, some 
truly remarkable realizations of single-mode non-Gaussian states have been recently carried out through 
the use of parametric amplification plus postselection lfT3l [T4l [TSl . Recently, by a generalization of the 
experimental setup used in Ref. ifTSll to a two-mode configuration, a method has been proposed ||3l for the 
generation of a certain class of two-mode photon-subtracted states. 

Here, in an analogous manner to that of ref. ||T3l , we propose a possible experimental setup for the 
generation of the degaussified states (eqs. ( 13.8b and ( 13.7b ). and of the squeezed Bell-like states (eq. (13.131) ). 
The scheme, based on a configuration of cascaded crystals, is depicted in fig. ( 13.5b . 

In the first stage, by means of a three-wave mixer, functioning as a parametric amplifier, a two-mode 
squeezed state \C)i2 = 5'i2(C)|07 0)i2 is produced. In the second stage, a four-wave mixing process takes 
place in a crystal with third order nonlinear susceptibility x'^'. We consider two possible multiphoton 
interactions, in the travelling wave configuration, described by the following Hamiltonian operators; 



r. 



states 




-t , * - - - 
KB^ia^a^ + Kgaia2a3 , 



(3.27) 



(3.28) 



where di (i — 1,2,3) denotes three quantized modes of the radiation field. The complex parameters ka 
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Figure 3.12: Scheme for the generation of the photon- subtracted squeezed state and of the photon-added 
squeezed state. Two nonlinear crystals are used in a cascaded configuration. The first x^^^ -crystal is part of 
a three- wave mixer, acting as a parametric amplifier for the production of a two-mode squeezed state. The 
squeezed state seeds the successive nonlinear process, a four-wave mixing interaction occurring in a x*^^'- 
crystal. A final conditional measurement reduces the multiphoton state to a photon-subtracted (or added) 
squeezed state | 5* )out. 

and K B are proportional to the third order nonlinearity and to the amplitude of an intense coherent pump field, 
treated classically in the regime of parametric approximation. The two-mode squeezed state seeds modes 1 
and 2; mode 3 is initially in the vacuum state | )3; mode 4 is the classical pump. Energy conservation and 
phase matching are assumed throughout. Let us remark that, due to the typical orders of magnitudes of the 
third order susceptibilities, the parametric gains are very small | | , | | ^ 1. The propagation (time 
evolution) in the crystal yields | ^I'j^'' ) = cxp{— it_ffj^''}| C)i2 |0)3 {L ~ A,B). Truncating the series 
expansion of the evolution operator at the first order in = — it kl^nq get 

l*^"")) «{1 + KAaia24}IC)i2|0)3 (3.29) 
«{1 + ^saI44}IC)i2|0)3 (3.30) 

Finally, a conditional measurement is performed on mode 3, consisting in a single-photon detection; a 
projection onto the state |1)3. The postselection reduces the states of eq. ( |3.29l l and eq. (13.30b to the states 
of eq. ( 13.7b and eq. (13.8b . respectively. It is worth noting that the low values of the parametric gains do not 
affect the implementation of the process. In fact, it is analogous to require low reflectivity of a beam splitter 
to generate photon-subtraction (addition) by using linear optics. 

Regarding the production of the squeezed Fock state (eq. (13.1b ). it can be generated, in principle, by 
seeding a parametric amplifier with a single-photon states in the two modes. 

Let us now turn to the experimental generation of the squeezed Bell-like states (eq. ( 13.13b ). They can 
be engineered by using the same setup illustrated in fig. (13.5b . and by simultaneously realizing inside the 
nonlinear crystal the processes corresponding to the interactions in eq. (13. 27b and eq. (13. 28b . In this case 
the fundamental requirements are that of energy conservation and phase-matching for each multiphoton 
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interaction must hold simultaneously at each stage. This condition can be satisfied by suitably exploiting the 
phenomenon of birefringence in a negative uniaxial crystal 1821 . In particular, the following set of equations 
must hold: 

K1'=' ^ k°"^ + k°"' + kl""* (3.31) 

where ujj and fc^ (j = 1, 2, 3) represent the frequencies and the wave vectors of the quantized modes with 
polarization A; ^Ij and {j ~ 1, 2) represent the frequencies and the wave vectors of the classical pump 
fields; the superscript ord and ext denote, respectively, the ordinary and extraordinary polarizations for the 
propagating waves. A coUinear configuration is assumed for the geometry of propagation inside the crystal. 
Then, at fixed uji and uj2, the energy conservation relations, the type-II phase matching condition in eq. ( 13.311 ). 
and the type-I phase matching condition in eq. ( |3.32t can be, in principle, satisfied by a suitable choice of 
UJ3, rii, fl2, and of the phase-matching angle between the direction of propagation and the optical axis. 
Various examples of such simultaneous multiphoton processes have been demonstrated both theoretically 
and experimentally lU [83] [84l US] . The final conditional measurement on mode 3 yields the superposition 
state 

\^i) « KAai 02^12(0 |0, 0)i2 + iiBa\alSi2iO\0,0)i2 (3.33) 

By applying a standard Bogoliubov transformation and after a little algebra, it is straightforward to show 
that the superposition state (eq. (13.33b ) reduces to the squeezed Bell-hke state (eq. (13.131) ) if 

ci = - (e-*'^KB tanh(r) + e''^ Ha) 

C2 =kb + e^^'^KA tanh(r) (3.34) 

The latter conditions can be successfully implemented by observing that the complex parameters ka and 
Kb can be controlled to a very high degree by means of the amplitudes of the external classical pumps. 
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Chapter 4 



Teleportation with "truncated" 
Gaussians and squeezed cat-like 
resources 

We introduce, in section \4~T\ the class of squeezed symmetric superpositions of Fock states ifJTl as a further 
generahzation of the squeezed Bell-like states (see eq. ( I3.13l l). and apply to them the optimization procedure 
defined in ref. Il24ll and used in chapter [3j obtaining a maximal fidelity of teleportation of a given input, 
over the superposition parameters defining the character of the resource. We find that the optimal resources 
for teleportation in this class of states are necessarily constrained to be second-order "truncations" of the 
two-mode squeezed Gaussian states (see eq. (I3.12l i). 

In section |4~2l we introduce the class of squeezed cat-like states; optimize this class of states for maximal 
fidelity of teleportation over the available parameters, and compare results with those obtained in the previous 
section. 

4.1 Truncated squeezed Gaussians: Symmetric Superpositions of Fock 
States and Bell-like states 

In chapter |3] the squeezed Bell-like states have been exploited as non-Gaussian entangled resources that 
generalize the search for an optimal teleportation fidelity among various classes of input states, which these 
states interpolate. The squeezed Bell-like states (eq. ( I3.13l l) can be formulated as the application of the two- 
mode squeezing operator (eq. (|1.41t ) to a general superposition of the two-mode vacuum and two-mode 
single-photon Fock state (having two photons one in each mode), similar to a Bell state of discrete variables 
(eq. CH). 

The optimization of teleportation fidelity over the superposition parameters yields a remarkable enhance- 
ment in the success probability of teleportation for various input states 1241 using squeezed Bell-like states. 
In this section, we produce a further generalization on the squeezed Bell-like states deemed the Squeezed 
Symmetric Superposition of Fock states 1371; and show that the optimal fidelity choice for this new class 
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of states (and all the squeezed Bell-like states) can be regarded as truncations (on the photon number) of 
two-mode squeezed vacuums. 

4.1.1 Definition: truncated Gaussians 

The squeezed Bell-like states, having only two superposition components and the constraint of normalization, 
can be parameterized in alternative ways. A parameterization based on the first order {n = 1) truncation (see 
eqs. ( |1.43t and ( |3.12| i) of the two-mode squeezed vacuum state is possible. Take a trivial generalization of 
the aforementioned expressions, 

I - (r + s) )ab = SAB{-r) Sab{~s) \0,0) ab 

OO 

= SAB{~r) (cosh(s))-i (tanh(-s))" \n,n)AB (4.1) 

n=0 

and delete all the terms of order n > 1. The normalization of the state holds for all s; genuine squeezing 
is performed with parameter — r, however. The only additional complication with respect to the former pa- 
rameterization is the relation between both; equations involving trigonometric functions of 5 and exponential 
functions of s. Thus, the squeezed Bell-like state can necessarily be regarded as the first-order truncation of 
some two-mode squeezed vacuum. 

A further and obvious generaUzation of the squeezed Bell-like states lies in considering squeezed, sym- 
metrical Fock states of higher number in the superposition. Let us consider the squeezed superposition 

|*)S5F= Sab{0 {cl + c\ + cly"^ 

X {co|0, 0)as + e'^^cill, 1)ab + e'"- C2\2 , 2) ab] (4.2) 

where c, and 9i are real constants and phases, respectively. This class of states, termed the squeezed 
symmetric superposition of Fock states ||J7) is not necessarily a second-order {n = 2) truncation of eq. (14.1b : 
given that there are three components of the superposition and only one normalization constraint. The second- 
order truncations of two-mode squeezed vacuums are special cases of the squeezed symmetric superposition 
of Fock states (SSSF). However, the choice of superposition parameters for the SSSF state that interests us 
most corresponds to an optimal fidelity of teleportation for a given input state. 

A convenient parameterization of the Ci coefficients in eq. ( 14.21 ) is provided by the hyper-spherical co- 
ordinates in three dimensions: cq = cos((5i), ci = sin((5i) cos((52), C2 = sin((5i) sin((52). It is expected 
that having an additional free superposition parameter will allow us to find an optimal teleportation fidelity 
greater than the optimal fidelity obtained using the squeezed Bell-like state; for the latter is the special case 
of the SSSF, for (52 =0. 

The two-mode characteristic (eq. ( |1.54| ))function xssf of the SSSF state ( |4.2| i can be calculated easily. 
Let us recall eqs. ( 11.551 ) for the pure state characteristic function; eq. (11.46b for the Bogoliubov transfor- 
mation effected by two-mode squeezing on a displacement operator; eq. ( 13.3b for the matrix element of the 
displacement operator. Following a procedure similar to the calculation of the characteristic function for 
the degaussified resource states (and for the squeezed Bell-like state, see subsection 13. 1.2b we obtain, for 
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?i = ^2 = 0; 

g-i/2(i5:4i^+ic^i^) 



[c^ + c( + ci) 

2 2 I 

X EE <^™Tr(^A^Br""^ir™'(ieAi')^^r'"^(iCBi') (4.3) 



where and are related to and by the Bogoliubov transformation for the displacement operator 
arguments in eq. ( 11.46b : and Lm"™^ is the associated Laguerre polynomial. 



4.1.2 Optimized teleportation fidelities with truncated Gaussians 

We proceed now to calculate the teleportation fidelities for ideal CV teleportation (see eq. (12.551 1) for coherent 
state and Fock state inputs, using the non-Gaussian resource states outlined above as resources (for the 
description and characteristic functions of the input states see subsection (13.1.4b . 

For the SSSF state, and the squeezed Bell-like state, we have calculated the teleportation fidelities and per- 
formed the optimization of said fidelity over the superposition parameters regarding each state. The squeezed 
BeU-like state is optimized over its sole free parameter 5 (see subsection 13.3.2b . The SSSF teleportation 
fidelity is in principle a function of parameters r, (f), 5i, 6i, 62, 6*2; but r is fixed as a technical capability 
indicator and the phase = tt is fixed as required by the conventions in the CV protocol of chapter |2] 
Fidelity is thus to be optimized over the set of superposition parameters V ~ {Si, 6*1, (52, 6*2}; 

Topt = max T{r , TT , Si , 9i , S2 , 02) (4.4) 

for a given class of input states. The optimization is performed numerically on the analytically calculated 
fidelities (reported in appendix lAl). The results of optimization are independent of the phases 9i and 6*2, 
which can be set to zero in consequence. In this manner the SSSF resource is sculpted for the teleportation 
task at hand. 

The procedure leads to a very remarkable result: Optimization over the set V yields maximal parameters 
for an SSSF state that is a second-order truncation of the two-mode squeezed vacuum state (that of eq. ( 14.11 )). 
The optimal resources for CV teleportation form a subset of the class SSSF having the form 

W)ssF - [1 + (tanh(,s))2 + (tanh(s))4]-i/2 

X SAB{-r) {|0,0)Ai3 + tanh(s)|l, + {i&n\i{s)f \2 , 2) ab] (4.5) 

where, in accordance with the optimization performed in eq. ( 14.4b . the real parameter s is fixed to the optimal 
value s ~ s. We have performed a change of parameterization that simplifies the description of the optimal 
resource: nevertheless J-opt{s) = maxg J^{r, s). 

In fig. 14.1.21 we report the behavior of Topt as a function of r for coherent state inputs (Panel I) and the 
single-photon Fock state input (Panel II) with a SSSF state as the entangled resource. For comparison, 
the curve for optimal fidelity corresponding to a squeezed Bell-like resource and the curve for fidelity 
corresponding to a two-mode squeezed vacuum are shown. The parameterization chosen for the squeezed 
Bell-like resource is that of the Gaussian truncation illustrated before (see eq. ( 14.1b and related discussion): 
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thus, an optimal "superposition" squeezing s is given for both non-Gaussian resources. We observe that, as 
expected, a further enhancement of the fidelity is obtained for the SSSF resource compared to the squeezed 
Bell-like resource. Even if the number of "free" parameters (one, s) appears to be the same for both non- 
Gaussian resources, the extra dimension in the Hilbert space spanned by the superposition allows an increase 
in fidelity. 




Figure 4.1: Optimal fidelity of teleportation J-opt, as a function of the squeezing parameter r, with = tt 
for truncated Gaussian resources. In Panel I, we report the fidelity for the teleportation of input coherent 
states I /3 ) using different two-mode entangled resources: ( a) Gaussian two-mode squeezed state (full line); 
fb) squeezed Bell-like state (dashed line); (c) squeezed symmetric superposition of Fock states (long-dashed 
line). In Panel II, we report the fidelity for the teleportation of input single-photon Fock state 1 1 ) using 
different two-mode entangled resources: ( a) Gaussian two-mode squeezed state (full line); (b) squeezed Bell- 
Uke state (dot-dashed line); (c) squeezed symmetric superposition of Fock states (double-dot-dashed line). In 
plot I the value of /3 is arbitrary. The insets in both Panels give the maximal values of the parameter s ~ s 
as a function of r for a given entangled resource and fixed input state. The plot styles are chosen as specified 
above. 

4.1.3 Entanglement, non-Gaussian character and Gaussian affinity for truncated 
Gaussians 

In section 13.41 to understand the properties of the optimized teleportation resources, we have investigated 
three quantities: the von Neumann entropy Ey^ (eq. ( 13.241 )) as a measure of the amount of entanglement in 
the resource; the non-Gaussianity d^Q (eq. ( 13.251 )). to provide a measure of the non-Gaussian character of the 
resource; and the squeezed vacuum affinity Q (eq. ( 13.261 )), in order to determine the degree of resemblance of 
the resource to the closest two-mode squeezed vacuum. 

In fig. ( 14.1.31 ), we plot the entanglement measure for pure states E,jn (eq. ( |3.24| l); for the optimized non- 
Gaussian entangled resources above mentioned and used as teleportation resources, and for the two-mode 
squeezed vacuum. All the curves exhibit very similar behaviors. At fixed r, and for a given input state, the 
optimized SSSF state is the most entangled; the optimal resource for the teleportation of a single-photon Fock 
state having more entanglement than the resource for the optimal teleportation of a coherent state; this can 
lead to the conclusion that the non-Gaussian input requires more entanglement for optimal teleportation. 

The behavior of the von Neumann entropies across different classes of resources is in agreement with 
the behavior of the optimal fidelities. In fact, for a given input state, a non-Gaussian resource with a higher 
teleportation fidelity (in fig. ( 14.1.21 )) is associated with a higher amount of entanglement content for any r. 

It is not remarkable that the Gaussian truncation (see eq. ( 14.11 )) to higher order (ji = 2, optimal SSSF 
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state)) shows a higher von Neumann entropy that the Gaussian truncation to lower order (n = 1, optimal Bell- 
like resource). We can only speculate that a higher-order (n > 2) truncation will have a higher von Neumann 
entropy for the same squeezing r. What seems remarkable, is that for the same squeezing parameter r, the un- 
truncated Gaussian has the lowest von Neumann entropy. Before beginning to speculate as to the truncation 
order in which the von Neumann entropy begins to diminish, let us think again about what is meant by 
Gaussian "truncation". Inspecting eq. (14. It ; and on seeing the squeezing — (r + s) performed on the resource 
to be "truncated", we can only conclude that to be "fair" to the Gaussian resource, the optimal non-Gaussian 
resources used here (defined by their parameters — r and — ,5, see eq. i4.5\ ) should be compared to a two- 
mode squeezed vacuum of squeezing — (r + s) instead. We rephrase our claim: "truncating" the Gaussian 
resource Q raises the entanglement of the resource: from that associated to a squeezing r for a two-mode 
squeezed vacuum, to that associated to a squeezing r + s' for the same state. This added squeezing s' is a free 
superposition parameter to be chosen. From inspection of figs. (I4.1.2l i and (14.1.3b : for the optimal resources 
we have considered the added squeezing s (and its effect on entanglement) is of the order of magnitude of 
the squeezing r (and its associated levels of entanglement). 




Figure 4.2: Plot of the von Neumann entropy E^n as a function of r, with s = s (see the insets in 
fig. I4.1.2l i. corresponding to the following entangled resources: squeezed symmetric superposition of Fock 
states, squeezed Bell-like states, and two-mode squeezed vacuum states, optimized for the teleportation of 
coherent state inputs, and the single-photon Fock state input. The plot styles are chosen as specified in 

fig. (133. 

Next, we examine the behavior of the non-Gaussianity dnG, a measure of the difference between the 
resource and a reference Gaussian state with the same first and second-order moments. In fig. ( 14.1.3b . panel 
I: we plot dnG (eq. ( 13.251 )) for the SSSF resources and the squeezed Bell-like resources, optimized both for 
the coherent inputs and for the single-photon Fock input. The difference in non-Gaussianity between the 
coherent input and Fock input cases for r = is remarkable. While the resources are entangled even with no 
squeezing; it seems that optimization for very low r is just an attempt at (nearly) classical teleportation. The 



We do not necessarily think that "truncation" is a physically feasible operation on a Gaussian resource. 
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behavior of the pairs of dnc curves corresponding to the teleportation of the same input states (either coherent 
or Fock) for different resources (SSSF and Bell-like) is very similar; a pair of curves (the Fock state input's) 
crosses at a certain value of r. All the curves tend to the same small asymptotic value of dnc for very large 
r, as expected. Therefore, according to this measure of non-Gaussianity: the non-Gaussian resource with a 
higher teleportation fidelity (the SSSF state) does not exhibit a significantly higher non-Gaussian character 
for any r. 
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Figure 4.3: Plot of the non-Gaussianity dnC (panel I) and of the squeezed vacuum affinity Q (panel II) 
as a function of r, with s = s (see the insets in fig. ( 14.1.21 )). corresponding to the following entangled 
resources: squeezed symmetric superposition of Fock states, and squeezed Bell-like states, optimized for the 
teleportation of input coherent states and the single-photon Fock input state. The plot styles are chosen as 
specified in fig. ( |4.1.2| i. 

Lastly, we turn our attention to the Gaussian resource affinity Q (eq. ( |3.26l l). given as a measure of the 
maximum similarity of the resource to an ideal two-mode squeezed vacuum. In fig. (|4.1.3l l panel II: we plot Q 
for the SSSF and squeezed Bell-like resources, as always optimized for the teleportation of given input states. 
We can see a mirror of the behaviors of dnc for r ^ 0. All the curves have a high asymptotic value of Q for 
large r, a fact consistent with the results for non-Gaussianity dnc- Most importantly: For the same class of 
input state {either coherent or Fock), the optimized SSSF state possesses a greater squeezed vacuum affinity 
than that of the squeezed Bell-like state, for any r. This result is in complete agreement with the behaviors of 
the corresponding optimal teleportation fidelities and von Neumann entropies studied before. 

Comparing the behaviors associated with resources optimized for the efficient teleportation of the same 
input states, we conclude that the following hierarchy can be established: The enhancement of the teleporta- 
tion fidelity corresponds to an enhancement of the entanglement content, and to an enhancement of squeezed 
vacuum affinity. The measure of non-Gaussianity d^Q seems not to have a place in this hierarchy. This result 
shows that such a measure must be used and interpreted with attention, case by case. 

Further conclusions can be drawn from the above discussion: The squeezed truncated Gaussian resource, 
with truncation at n — 2 (SSSF state, with a four-photon term) is an entangled non-Gaussian resource 
suitable for sculpturing and optimization as regards CV quantum teleportation. We argue that this result can 
be generalized to higher orders (n > 2) in expansions of the two-mode squeezed vacuum. 

Therefore, we conjecture that in general, efficient entangled non-Gaussian resources for CV quantum 
teleportation should be characterized by a suitable balance between the entanglement content and the degree 
of affinity to squeezed vacuum states. 
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4.2 Squeezed cat-like states 

In this section, we define and use as teleportation resources a special class of two-mode squeezed superpo- 
sition. This time, instead of using (solely) Fock states for the superposition we use coherent states. Namely 
the vacuum (which is a coherent state) and an arbitrary separable coherent state with the same displacement 
on both it's modes. These can be thought of as more feasible and more amenable to manipulation than the 
superpositions made of the highly nonclassical Fock states. We remark that the engineering of states of such 
a form is made possible, in principle, by the recent successful experimental realization of their single-mode 
equivalents ifTSl . 

4.2.1 Definition: squeezed cat-like states 

The squeezed cat-like state is given by the two-mode squeezed superposition of coherent states 

\i^)sc = ^fSABiC) {cos(J)|0, 0)ab + e'*' sin(<S)|7, 7)as} (4.6) 

where the normalization factor is JV = il + e '''' sin(2 S) cos{9) j . Without the two-mode squeez- 
ing, eq. l4.6l becomes an entangled Schrodinger Cat state. The one and two-mode cat states have been proposed 
as qubits and as entangled resources (respectively) for quantum information processing; together with an 
universal set of operations for this purpose ||30| . 

The state in eq. ( 14.61 ) can be regarded as the substitution of a coherent state | 7 , 7 ) ab for the two-photon 
I 1 , 1 )ab, in the squeezed Bell-like state of eq. ( 13.13b . Thus producing a more classical two-mode squeezed 
superposition that is non-Gaussian and entangled even for zero squeezing r. 

Following the same reasoning used to calculate the characteristic function of the SSSF state and using 
mostly the same identities, plus the composition law for displacement operators (eq. dl.lSb ). the characteristic 
function for squeezed cat-like resource of eq. ( 14.6b can be easily calculated. For = it is given by 

+ (sin(J))2e2'i'"[''*(«A+?B)]) (4.7) 

where and are related to and £,b by the Bogoliubov transformation for the displacement operator 
arguments in eq. (11.46b . 

4.2.2 Optimized teleportation fidelities with cat-like states 

We can compute analytically the fidelity of teleportation using squeezed cat-like states as entangled resources, 
limiting the analysis, for simplicity's sake, to the teleportation of coherent state inputs (eq. (13.15b ). The 
maximization procedure is discussed in appendix |A] including the analytical expression for the fidelity. At 
fixed squeezing r, the fidelity in eq. ( IA.5b may be optimized over the real amplitude | 7 | (the only free 
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parameter remaining). Therefore; 



opt{r) = max J^sc {r, hi ) 

I 7 I 



(4.8) 



In fig. (14.2. 2l l. we show the behavior of Topt as a function of r for coherent state inputs, using an optimized 
squeezed cat-like state as a teleportation resource. For comparison, the optimal fidelities corresponding to 
the teleportation with Gaussian squeezed state resources and with optimized, squeezed Bell-like states are 
plotted. The optimized squeezed cat-like resources, for all r, yield a significant improvement of the fidelity 
with respect to the two-mode squeezed vacuum, but are less efficient than squeezed Bell-like resources; all 
three curves show the same behavior, including the same asymptotic behavior. The high value of 7 for 
r = is a result of the optimization trying to increase the usable entanglement that is present in a two-mode 
Schrodinger cat state. For much larger squeezing, the optimal amplitude 7 tends to a smaller, asymptotic 
value; this behavior is consistent with the behavior of the resource, which eventually becomes an EPR state, 
for r —!■ 00. 



opt 




Figure 4.4: Optimal fidelity of teleportation J^opt, as a function of the squeezing parameter r, with 4> = tt, 
using squeezed cat-like resources. The fidelity corresponds to the teleportation of input coherent states with 
a squeezed cat-like state as entangled resource (dotted line). For comparison, we plot the fidelities obtained 
with a squeezed vacuum resource (full line) and with a squeezed Bell-like resource (dashed line). The inset 
gives the optimal value of the parameter | 7 | = | 7 | as a function of r. 



4.2.3 Entanglement, non-Gaussian character and Gaussian affinity for cat-like states 

The behavior of the teleportation fidelity for the optimized, squeezed cat-like resource in fig. (14.2.2b is in 
agreement with the behaviors of its von Neumann entropy (eq. ( 13.24b ). non-Gaussianity (eq. ( 13.25b ). and 
squeezed vacuum affinity (eq. (13.26b ). We report these in fig. ( 14.2.3b and fig. (14.2.3b . panel I and panel II, 
respectively. 

The optimized squeezed cat-like resource is less entangled, less non-Gaussian (according to the dno 
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measure), and less affine to the squeezed vacuum than the optimized squeezed Bell-like state. The entangle- 
ment as measured by the von Neumann entropy is higher for the cat-like resource than for a Gaussian, but 
lower than that for a Bell-like state for all values of r. The optimized squeezed cat-like resource exhibits 
a markedly lower value of Q in comparison with the truncated Gaussian resources. This can be seen by 
comparing fig. ( 14.2.31 ). panel II with the corresponding plots shown in fig. ( I4.1.3I I, panel II; in fact, Q seems 
to have an asymptote at value 71% (approached from below). 
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Figure 4.5: Behavior of the von Neumann entropy E^jq as a function of r, with | 7 | = 7 (see the inset 
in fig. (I4.2.2I 1) for the optimized squeezed cat-like state (dotted line). For comparison, we also plot the 
quantities E^j^ corresponding to a two-mode squeezed vacuum (full line), and to a squeezed Bell-like state 
(dashed line). 




Figure 4.6: Behavior of the non-Gaussianity dnc (panel I) and of the squeezed vacuum affinity Q (panel II) 
as a function of r, with | 7 | = 7 (see the inset in fig. (14.2.21 )) for the optimized squeezed cat-like state (dotted 
line). For comparison, in panel I we plot dnQ for the optimized squeezed Bell-like state (dashed line). 
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Chapter 5 

Teleportation with noisy non- Gaussian 
resources 

In chapter[3] we defined the squeezed Bell-like states (eq. ( 13.13b ); a class of states that include all the other 
non-Gaussian resources studied in that chapter, as well as the Gaussian two-mode squeezed vacuum. We 
optimized the Squeezed Bell-like state for the teleportation of some inputs and compared them for their 
properties of entanglement, non-Gaussianity and Gaussian resource affinity with the other resource states 
studied in that chapter. In chapterlH we introduced the squeezed cat-like states iJTl (eq. ( 14.6b ) by substituting 
a coherent state term for the Fock state term in the superposition making up a Bell-like state, with the intention 
of creating a more feasible ifTSl . more resilient |fT9l (to decoherence) resource. 

In this chapter, we extend the analysis performed in chapters [3] and |4] to the realistic case of optimized, 
squeezed Bell-like and squeezed cat-like resources prepared or propagated in the presence of thermal noise 
(see section lZ4l i. for the teleportation of coherent state inputs ll36l[37l . For comparison purposes we perform 
the same analysis for Gaussian two-mode squeezed vacuum states, also prepared in the presence of thermal 
noise. 

We analyze the behavior of inseparabihty parameters developed for general CV contexts BTl l42l . for 
BeU-like and Gaussian resource states. For a practical measure of the disappearance (or appearance) of 
entanglement of teleportation, we consider the level of noise at which the resource passes the classical 
teleportation threshold 1391 l40ll of 1/2, for a fixed squeezing. This last analysis is performed for both 
the squeezed Bell-like resource and the squeezed cat-like resource, with the squeezed vacuum put in for 
comparison purposes. 

5.1 Fidelity of teleportation with mixed, noisy non-Gaussian resources 

The presence of thermal noise in the preparation of the resource states can be modelled by superimposing a 
pure resource state on thermal states (see section 12741 ) in modes A and B of the resource. The characteristic 
function of such a preparation (see eq. ( 12.411 )) is given by 

xfi^iU^B) - e-"---l«-l'-"-'-l«-l'xAB(^A,?B) (5.1) 
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where nth^A, nth,B are the mean photon numbers; thermal parameters associated with the modes A and B. 

The state of eq. ( 15. lb is mixed for nth,A,nth_B 0. xab is the characteristic function for the pure 
resource state thus superimposed over thermal states. We substitute in eq. ( 15.1b the characteristic functions 
of either the squeezed Bell-like state (eq. ( 13.14b ): the squeezed cat-like resource (eq. ( 14.7b ) or the two-mode 
squeezed vacuum (eq. (11.58b ) to produce the respective "noisy" teleportation resources. 

Substituting the resource of eq. ( 15.11 ) in eq. ( 12.261 ) with a coherent state input (eq. ( 13.15b ) and ~ 
gp = 1; we obtain the output for CV teleportation with a noisy resource. From the output state, it is 
straightforward to calculate the analytic expression for teleportation fidelity as a function of the adequate 
superposition parameters: 5 for the squeezed Bell-like state and | 7 | for the squeezed cat-like states. The 
squeezing parameter r is fixed and = tt in keeping with convention. Likewise are fixed the thermal 
parameters nth, a and nth,B- The analytic expressions of the fidelities for the noisy non-Gaussian resources 
in the teleportation of coherent states; squeezed Bell-like {Tgl^^ ( r , nth, a , nth,B , <5 )) and squeezed cat- 
like (J^sc' I nth, A , nth,B , I 7 I )) and an initial analysis of the optimization procedure performed on the 
aforementioned fidelities are reported in appendixlAl 

For example, for the mixed squeezed Bell-like state teleporting a coherent state, the fidelity is optimized 
over S, leaving other parameters fixed; 

^opt{r , nth, A, nth, a) = max J^j^^(r, nth,A , nth,B , 5) , (5.2) 

o 

The optimization yields an optimal angle Smax' that has a form not too dissimilar to that for ideal 
teleportation (eq. ( 13.21b ): 

4^'"'^ = \ arctan ( 1 + , ^ ) (5.3) 

2 \ I + nth.A + nth,B J 

which reduces to the pure state case if the thermal parameters nth, a = nth,B — 0. When this angle 
is substituted in the expression for the fidelity (eq. (lA.lb ) gives the optimal fidehty for a mixed Bell-like 
resource. 

In the case of mixed squeezed cat-like states, at given thermal numbers nth.A, nth,B and fixed squeezing 
parameter r, the optimal fidelity is defined as 

^optir, nth,A, nth,B) = max T^*^\r , nth,A , nth,BAl\) (5-4) 

1 7 

where the maximization is performed numerically for fixed values of r and nth = nth, a = nth,B- 

In fig. (15.1b : we plot the optimal fidelities for the mixed non-Gaussian resources as a function of squeezing 
for various choices of the thermal parameters nth.A = nth.B = nth- Along, we plot the fidelities for 
a Gaussian two-mode squeezed vacuum, used as a benchmark. We observe that, as expected, the fidelity 
decreases for increasing nth- the thermal noise sensibly reduces the success probability of teleportation, as 
it reduces the entanglement content of the resources. We can also observe that for similar noise nth the 
asymptote for large r is the same for all resources: the asymptote indicates the convergence of all resources 
to a mixed two-mode squeezed vacuum at high squeezing, with the "mixedness" putting an upper bound on 
teleportation fidelity even for this nearly ideal resource. 

An important observation is that in our plot, for the chosen realistic values for 7ith', the fidelity associated 
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with both non-Gaussian mixed resources never drops below the threshold of classical teleportation 1391 with 
maximal fidelity T'^^^ = 1/2. The ability of the resource to keep the fidelity above this benchmark value of 
1/2 will be made a practical measure of resilience in presence of noise; in the next section. 

At fixed squeezing and thermal parameters, the non-Gaussian resources always perform better than 
the Gaussian two-mode squeezed vacuum. Furthermore, mixed squeezed Bell-like states have a higher 
performance than the mixed squeezed cat-like states. We remark that the better performance of the mixed 
Bell-like states over the mixed cat-like states has been demonstrated for this simple analysis of environmental 
noise. Other sources of noise (and decoherence) have not been studied. 




Figure 5.1: Optimal fidelity J-qja, associated with the teleportation of input coherent states with mixed 
squeezed Bell-like resources (panel I) and mixed squeezed cat-like resources (panel II), as a function of 
the squeezing r, for several choices of the thermal parameters nth, A = nth,B = nth- The curves representing 
the fidelities associated with mixed non-Gaussian entangled resources are plotted with the following plot 
style: nth = (full black line), nth = 0.05 (long-dashed line), nth = 0.10 (double-dot dashed line), and 
nth = 0.15 (dotted line). For comparison, in both panels, we also plot the fidelities associated to the mixed 
squeezed vacuum entangled resources, with rith = (full gray line), nth — 0.05 (dashed line), nth = 0.10 
(dot-dashed line), and nth = 0.15 (long-dotted line). The inset in panel II gives the maximal value of the 
parameter | 7 | = 7 as a function of r, with nth = 0, 0.05, 0.1, 0.15 (same plot style as for the fidelities). 

We have chosen the mixed Gaussian state to be a reference resource (see subsection |3.3.2| i for the 
evaluation of the performance of the squeezed Bell-like state in teleportation. To develop this relationship 
between resources further; we define the relative fidelity (see eq. ( 13.23b ) given for mixed squeezed Bell-like 
state and referring to the mixed Gaussian state 

^T-(c). ^ ^''''Kr,nth,nth,Stax'') - J^^''\r,nth,nth,0) 

^Kptir.nth) = TT. (5.5) 

J^(''>{r,nth,nth,0) 

where T^'^^{r, nth, nth, 0) is the fidelity for the mixed Gaussian resource and JF(^'(r, nth, nth, Smax^) is the 
fidelity for the optimized noisy Bell-like resource (eq. dA.U ). 



In fig. (15.1b we show the behavior of ATjj'^^^{r,nth) as a function of r. We see that the percent gain 
in the fidelity for nth = essentially coincides with the one defined in the absence of thermal noise 
(fig. (13.3.2b ). It is also evident that for ever increasing nth, both the resources will converge asymptotically 
to Aj^Qpj(r, nth) = 0. This signals a similar and very poor performance in teleportation for both resources; 
corresponding to a separable, mixed, thermal state of two-modes. 

We have shown that the mixed non-Gaussian resources; although with reduced fidelity with respect to 
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Figure 5.2: Relative fidelity for mixed Bell-like resources with respect to mixed Gaussian resources for 
coherent state inputs; AT^'^f{r,nth) as a function of the squeezing r, for several choices of the thermal 
parameter nth, a = nth g = "-t/i- "-i/i = (full black line), nth = 0.05 (long-dashed line), nth = 0.10 
(double-dot dashed fine), nth = 0.15 (dotted line). 



the ideal resources, are better choices than the mixed Gaussian two-mode squeezed resource. Both the non- 
Gaussian resources studied here also provide acceptable (always better than the Gaussian) values of the 
fidelity for technically feasible squeezing r and realistic values of thermal mean photon number nth- 



5.2 Inseparability criteria and classical teleportation: Bell-like, cat- 
like and Gaussian resources 

In ref. Il24l it has been shown that the entanglement, that is a fundamental requirement for the efficient 
implementation of nonclassical teleportation protocols, is remarkably enhanced in the Bell-like state when 
compared to the Gaussian squeezed vacuum. Even though the pure squeezed Bell-like state has entanglement 
for r = 0; the introduction of noise (and of a mixed nature) induces a lowering of the amount of entanglement 
in the Bell-like state. It is then necessary to check that, for fixed values of the thermal parameter nth, a 
sufficient amount of entanglement still survives and is useful for teleportation, above a classical threshold 
value. 

In order to avoid the not-yet-solved problem of computing the amount of entanglement by means of a 
measure appropriate for all mixed states; we can exploit an inseparability criterion, based on the condition 
of positivity under partial transposition (PPT criterion) lf86l [87l [88l 1411 l42l . For the purpose of knowing 
where (for which values of r and nth) a state becomes inseparable and possibly useful for quantum telepor- 
tation, we assume this criterion to be just as useful as a proper entanglement measure. 

In our case, the PPT criterion can be expressed through an inequality for a inseparability parameter A 
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involving only second order statistical moments: 



A = {a\aA) {a^B^B 



aA as , 



< 



(5.6) 



which are straightforward to calculate, as the characteristic function of a quantum state is the generating 
function for the statistical moments of the state (see eg. (11. 571 )). 

Let us recall that the inequality in eq. ( 15.61 ) is a sufficient inseparability condition for non-Gaussian states, 
being necessary and sufficient for Gaussian states. In fig. (15.21 ): we plot the behavior of A as a function of r 
for several choices of the thermal parameters nth. a ~ nth,B ~ ^th, both for the mixed squeezed Bell-like 
state and for the mixed squeezed vacuum state. 
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Figure 5.3: Behavior of the inseparability parameter A as a function of the squeezing r, for several choices of 
the thermal parameters nth,A = nth,B = nth- Plots are for the mixed squeezed Bell-like state, with i^th — 
(full black line), nth = 0.05 (long-dashed line), nth = 0.10 (double-dot dashed line), nth = 0.15 (dotted 
line); and for the noisy squeezed vacuum state, with nth = (full gray line), nth = 0.05 (dashed Une), 
nth = 0.10 (dot-dashed line), nth = 0.15 (long-dotted line). 



Recall again that the pure Bell-like state is already entangled for r = and 6^0, 7r/2 (see fig. ( 13.4.11 )). 
For realistic values of nth (0 to 0.15) we see the mixed Bell-like state exhibiting entanglement. On the other 
hand, for nth > 0, the mixed squeezed vacuum state has A > at sufficiently low values of r; i.e. it 
becomes separable. Specifically, for the mixed Gaussian state the threshold value for separability is 

Though we have defined no measure for entanglement, and in particular, no measure of entanglement 
useful for CV teleportation; we have in the two-mode squeezed vacuum a practical reference resoMrce, with 
a "maximum classical fidelity" ^ = 0.5 that is widely accepted as a practical threshold 1391 l40l for CV 
teleportation of coherent states. We can assume for practical purposes that quantum teleportation becomes 
impractical when thermal noise causes fidelity to go under this value. 

Therefore, for a resource that teleports coherent states, we define the classical threshold value of the 
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thermal parameter n['f^'^^ (r): At a fixed r and for nth = "-i^'"'' (j')- the optimal fidelity Topt = ^S's^ = 1/2. 
The optimal fidelity goes above the 1/2 threshold in a smooth manner for nth < nfh^'^ (r), and goes smoothly 
below the threshold 1/2 for nth > "th*^('')- For this it is assumed (and found for the resources used here) 
that the optimal fidelity is an analytical function of r and nth- 

In fig. 15.21 we plot nl^'"'' as a function of r for the optimized, mixed squeezed Bell-like resource, the 
optimized mixed squeezed cat-like resource and the mixed squeezed Gaussian resource, shown as a reference. 
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Figure 5.4: Behavior of the classical threshold value n^^ of the thermal parameter as a function of r, for 
mixed squeezed Bell-like resources (dashed line), mixed squeezed cat-like resources (dotted line), and mixed 
squeezed vacuum, shown as a reference (full line). 

We see that the threshold value associated to non-Gaussian resources is larger than that associated 
to Gaussian ones, and in any case sensibly larger than the reasonable realistic values (i%th < 0.15) we 
have considered. Even for r = the non-Gaussian resources, always under realistic conditions, show 
entanglement allowing teleportation fidelities greater than the classical threshold. On the contrary, the fidelity 
associated with the mixed squeezed vacuum state, at low values of r falls below the classical threshold, 
unless optimization is performed on the local degrees of squeezing, without modifying the Braunstein- Kimble 
protocol 

ids) 

It is remarkable that the squeezed Bell-like state has a consistently higher value of nj,^ ' than the squeezed 
cat-like state. Taking this result at face value, it can be said that the Bell-like state is more resilient in noisy 
environments than the squeezed cat-like state, at equal r. Note lastly that for mixed squeezed vacuums, 
nflf^ (r) coincides with the threshold value for separability, that is nf^^^ (r) = ni^"^^ ir) = ■ 

Compare fig. ( 15.11 ) with fig. (|5.2| i for the squeezed Bell-like resource. In the former, increasing nth 
makes for a marked decrease in Ajr^pj(r, nj/i) for all r. In the latter, the asymptotic limit of 'f^th'^ if), for 
both Gaussian and non-Gaussian resources is the same. Therefore, we have a value of nth ~ 0.5 beyond 
which no reasonable amount of two-mode squeezing r will avail to produce a useful resource for quantum 
teleportation out of either a Gaussian resource or a Bell-like resource of any kind. Fortunately, this value of 

'More useful than a separable pure Gaussian state in any case. 
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nth is considered huge for a realistic experimental setting. 



80 Chapter 5 Teleportation with noisy non-Gaussian resources 



Chapter 6 



Conclusions 



We have developed a Wigner's characteristic function based formalism that allows for the intuitive represen- 
tation of the formalism for CV quantum teleportation for any combination of resource and input states for 
which these characteristic functions exist. The formalism is not confined to the original protocol; furthermore, 
it is possible to introduce further and more complicated steps to the protocol; measurements different from 
homodyne detection of EPR states and the mixing in of external modes representing noise, photon addition 
and photon subtraction. We have further shown that an analysis of the expression for the teleportation fidelity 
makes for a trivial derivation of results for resource states that are simple mixtures of pure states. The 
characteristic function formalism, as it stands, can easily accommodate most of the conceivable operations 
and measurements involved in CV teleportation, together with simple modifications to the protocol. It is 
then, possible, to use characteristic functions to "model" other quantum information protocols in a convenient 
manner. 

We have presented a thorough comparison, with regard to the performance in continuous-variable quan- 
tum teleportation, between standard Gaussian, wholly non-Gaussian (two-mode squeezed Fock state) and 
degaussified resources (such as photon-added and photon-subtracted squeezed states) and a new type of 
sculptured resource that interpolates between these states and can be optimized because it depends on an 
extra, relative-phase, independent free parameter in addition to squeezing. These sculptured non-Gaussian 
resources we have named squeezed Bell-like states: They hybridize discrete single-photon pumping, coherent 
superposition of Bell two-qubit eigenstates, and CV squeezing; thus including the above mentioned resources 
as "special cases". The maximization of the teleportation fidelity (an analytical expression which is valid 
for a given input state) is made with respect to the free parameter. Therefore we have produced a fidelity 
that is optimal for the class of Bell-like states including, most importantly, all the "special case" resources. 
Understanding the enhancement yielded by sculptured squeezed Bell-like resources in teleportation success 
is possible when certain properties of the resources used are studied and compared jointly. The optimized 
squeezed Bell-like states are those states that are as close as possible to the simultaneous maximization of 
entanglement, non-Gaussianity, and affinity to the two-mode squeezed vacuum. 

We have proposed a method for the experimental generation of squeezed Bell-like states using a cascading 
setup of second and third-order non-linear crystals. Given the nature of the squeezed Bell-like states, this 
method is to be additionally regarded as an alternative method for the generation of degaussified resources. 

We have further generalized the Bell-like states by including an additional four-photon term in the 



81 



82 Chapter 6 Conclusions 



superposition, producing squeezed symmetric superpositions of Fock (SSSF) states. These resources, with 
an added dimension for sculpturing, can be utilized to an optimal fidelity of teleportation that is a notable 
improvement even over the optimal Bell-like state fideUty. We have considered the Bell-like states as first- 
order truncations of a r+s squeezing Gaussian state where r is genuine two-mode squeezing and s depends on 
the character of the Bell-like state; we found, remarkably, that the optimal SSSF resource is to be formulated 
as the second-order truncation of the same r + s squeezing Gaussian states. We have argued that a Gaussian 
"truncation" should be compared to the r + s squeezing Gaussian instead of the r squeezing Gaussian. The 
analysis of the properties for these states, the Bell-like and the SSSF shows that the optimal teleportation 
resource comes close to simultaneously maximizing entanglement and affinity to the two-mode squeezed 
vacuum. Not so for the non-Gaussianity. 

We have introduced, and optimized for teleportation a class of two-mode squeezed cat-like states; su- 
perpositions of coherent states that are then two-mode squeezed. These states have been optimized for 
teleportation over their coherent amplitudes, obtaining performances and entanglement values that are higher 
than those for Gaussian resources; but lower than those for Bell-like states. These entangled states are 
interesting because of their, in principle, greater feasibility ifTSl in comparison with Bell-like and SSSF 
states, which involve superpositions of few-photon Fock states. 

Given the above results, we state that the optimization of the ideal CV teleportation protocol with non- 
Gaussian resources necessitates only the formulation of more general and more complicated non-Gaussian 
resources. The only difficulty to be incurred in the generalization of non-Gaussian resources Ues is in the 
analytic calculation of the fidelities for the purpose of optimization of the same. 

Further optimization is in principle possible with respect to the local parts of the resource states, in 
analogy to the case of standard Gaussian resources ||89l . One could think of extending the sculpturing to the 
entire basis of Bell states, to generate entangled non-Gaussian resources that can never be reduced to proper 
truncations of Gaussian squeezed resources. Such fully sculptured resources might allow for the further 
enhancement of the teleportation success due to the presence of a larger number of experimentally adjustable 
free parameters in addition to squeezing. Fully sculptured states could be applied to hybrid schemes of 
teleportation combining continuous-variable inputs with discrete-variable resources and viceversa. In this 
framework, a particularly appealing line of research would be to look for modified schemes of teleportation 
beyond the standard Braunstein-Kimble protocol, to be realized by generalized measurements in combination 
with state-control enhancing unitary operations. 

We have studied the efficiency of CV teleportation of input coherent states using, as resources, squeezed 
Bell-like and squeezed cat-like states, prepared, superimposed over initial thermal states that represent ther- 
mal noise. Thus we have defined general non-Gaussian, realistically mixed entangled states. We have 
shown that, although the thermal noise strongly affects the success probability of teleportation, the resource 
provided by the optimized mixed extensions of the non-Gaussian squeezed Bell-like and squeezed cat-like 
states guarantee a sufficiently high fidelity, for realistic values of the average thermal photon number nth 
and of the squeezing r. A better performance is assured always, with respect to the mixed extension of the 
Gaussian twin beam. 

We have calculated some simple expressions for the characteristic function of the teleportation output for 
different models of nonideal Bell measurement. The analysis of "noisy" teleportation that we have performed 
can be generalized using the above mentioned results, while taking into account the decoherence induced by 
the propagation of the resource in noisy channels. 
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Finally, the present discussion could be extended to other types of quantum information tasks and pro- 
cesses besides teleportation. For instance, it would be interesting to investigate the comparative effects of non- 
Gaussian inputs and non-Gaussian resources in schemes for the generation of macroscopic and mesoscopic 
optomechanical entanglement ||90l. A further goal will be the study of optimal teleportation with non- 
Gaussian resources of two-mode and multimode states 19 1 1 . 



Appendix A 

Analytic expressions of Teleportation 
fidelities 



In this appendix we report the explicit, analytic expressions for the teleportation fideUties calculated in 
chapters [3] and |5] We have deemed these expressions both too cumbersome and unnecessary for the purposes 
of our exposition regarding the CV teleportation protocol with Non-Gaussian resources; relying instead on 
the plots accompanying the exposition for the graphical representation of said fidelities. 



A.l The fidelities: noisy squeezed Bell-like states with coherent state 
inputs 

The fidelity !F'^''\r , nth, a , nth,B , <5) (superposition phase 6' = 0; squeezing angle (f> — n) for the CV 
teleportation of a coherent state input of arbitrary amplitude, using as a resource the mixed squeezed Bell- 
like state (see eq. ( 15. H and eq. ( |3.14t reads 



Wc). „ „ .^ l + e'^fth + e^'-/,^ + e^^U cos(2^) + [1 + V,,] sin(2^) 
'[r, nth,A, nth,B , 0) = ^-2ru — , „2r--f 1.3 

where 



11 + e^^/t, 



fth = I + nth,A + nth.B (A.2) 

The fidelities for the states that are special cases of the Squeezed Bell-like state (see section ITTl i can be 
obtained by the substitution the appropriate values of the 5 parameter in eq. (lA.ll i and the appropriate levels of 
noise nth, a and nth,B- For the case in which 5 — Q and nth,A — nth,B=0, the fidelity equals the well-known 
result 

1 + e ^'^ 



holding for pure two-mode squeezed vacuum resources and coherent state inputs 113911 . 

In the absence of thermal fields, for nth, a = nth,B = 0; we have fth = 1, and the fidelity reducing to the 
pure resource fidelities studied in chapter[3] Notice that, in the limit of large thermal parameter, the fidelity 
of teleportation in eq. ( lA.lb is strongly suppressed and vanishes asymptotically. 
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A.2 The fidelities: noisy squeezed cat-like states with coherent state 
inputs 

We report the analytical expressions for the fidelities of teleportation of input coherent states using pure and 
mixed squeezed cat-like states as resources. By putting, as usual, the phases 6 — and 4> ~ tt, the fidelity 
for the squeezed cat-like state (eq. (I4.6l l) reads 

. , ^ . cos^((5) + e 1 + ^'^'- s\n'^{6) + e"!"*"! (ei + = ^" + eA + ''■"' ) sin((5) cos((5) . 
^^'^("'^'^^ = (l + e-^-)(l + e-l7Psin(2^)) ^^'^^ 



It is worth noticing that for 5^0 and/or 7^0, the fidelity in eq. ( IA.4b reduces to the well known 
expression of eq. (IA.3b . A preliminary numerical optimization procedure for eq. (I A. 4b allows us to fix the 
parameters arg7 and 5 to the values arg7 = and 5 = ^, leading to the simplification of the above fidelity 
to the expression 



^sc {r , I 7 I) = (i+'-';')(7+e-l7P) ^"^'^^ 

Thus, at fixed squeezing r the maximization can be carried out with the real amplitude | 7 | being the only 
free parameter We have Topt ~ max| ^ | Tsc {1" , I 7 I)- In the limit of zero squeezing (r = 0), the value of 
the optimal fidelity becomes jc-^pt = [4 (^2 - 1)]-^ ~ 0.6035, for | 7 | = ln^/^(V2 - l)"^ ^ 1.3276. 

Lastly; consider the "noisy", mixed squeezed cat-hke resource with fixed angle 5 = j and phases <f) = n, 
6* = 0, arg7 = 0. The fidelity J-g^'^ (r , n^^t 1 , n^^i 2, 1 7 I) for the mixed squeezed cat-like state is given by 

^SC , -^'.a , n*.,2, 1 7 I) = e--(l + e-/*.)(l + e-I.P) ^^'^^ 
where /f/i is given by eq. ( |A.2b . 
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